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ABSTRACrr 

This  thesis  consists  of  a  three-dimensional  numerical 
analysis  of  the  Jet  Aircraft  Hush  House  located  at  Naval  Air 
Station  Jacksonville,  Florida.  Utilizing  the  PHOENICS  Code 
allows  for  the  determination  of  the  aerothermal 
characteristics  including  velocity,  pressure,  enthalpy, 
turbulent  kinetic  energy  and  the  dissipation  rate  of  turbulent 
kinetic  energy  in  the  facility  during  testing  of  the  U.S. 
Navy's  F-4  (Phantom  II)  J-79-GE-8  gas  turbine  engine  with 
afterburner.  How  and  by  what  method  PHOENICS  arrives  at  this 
solution  is  discussed.  Of  greatest  importance  is  the  resulting 
"behavior"  of  the  aerothermal  system. 

Problems  encountered  using  the  PHOENICS  Code,  resulting 
numerical  solutions  to  the  particular  facility,  comparison  to 
actual  test  data  and  recommendations  for  further  applications 
of  the  PHOENICS  Code  are  presented. 
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I  -  INTRODUCTZC3N 


New  and  overhauled  aircraft  gas  turbine  engines  undergo 
stringent  performance  testing  before  returning  to  operational 
use.  To  allow  for  safe  and  effective  monitoring  of  the  engine 
and  aircraft,  gas  turbine  Test  Cell  and  Hush  House  facilities 
are  utilized  to  perform  monitoring  of  key  temperatures, 
turbine  speeds  and  flow  rates  associated  with  the  engine.  The 
test  facility  additionally  provides  distortion-free  air  flow 
to  the  engine  inlet  (eliminating  rain  and  wind  effects) ,  for 
lower  exhaust  temperatures  exiting  from  the  facility  and 
reduced  noise  pollution  to  the  surrounding  area.  Furthermore, 
these  facilities  allow  for  final  engine  "trim"  settings  to 
ensure  that  maximum  power  is  available. 

Currently  the  test  facilities  designed  and  built  for  the 
Navy  vary  widely.  The  various  facilities  include; 

1.  The  Run-up  pad — where  a  tied  down  aircraft  conducts 
power  checks  in  an  open  air  environment. 

2.  Sound  suppressor — located  at  the  Run-up  pad  to  reduce 
noise . 

3.  Open  Test  Stand  (OTS) — located  outdoors  that  supports  an 
aircraft  engine  and  has  portable  instrumentation 
attached . 

4.  Hush  House — a  totally  enclosed  facility  allowing  around 
the  clock  testing  of  aircraft. 

5.  Test  Cell  {TC)--an  enclosed  facility  allowing  testing  of 
various  engine  sizes  out  of  the  airframe. 


Design  engineers  have  used  experience  and  common  sense  to 
build  these  facilities  with  varying  degrees  of  success. 
Aerothermal  monitoring  supplements  the  design  of  the 
structures. 

With  the  tremendous  technological  advances  made  in  recent 
years,  much  higher  temperatures  and  thrusts  are  available 
compared  to  similar  size  and  weight  engines  of  just  a  few 
years  ago.  Considering  that  nearly  67%  of  the  Navy  Test  Cells 
and  Hush  Houses  are  at  least  25  years  old  it  is  evident  that 
a  new  generation  of  standardized  facilities  meeting  federal, 
state  and  local  noise,  pollution,  security  and  safety 
standards  is  required.  [Ref.l:p.33 

The  Navy  Environmental  Health  Center  cited  total  claims 
against  the  government  in  1981  at  $114  million  due  to  physical 
injury,  hearing,  or  respiratory  damage  to  military  and 
civilian  personnel.  Additionally,  the  1981  ruling  by  the 
federal  court  system  that  military  gas  turbine  test  facilities 
must  comply  with  Section  233  of  the  Clean  Air  Act  accelerated 
the  need  for  modern  test  facilities  meeting  the  above 
requireme:ts  plus  those  projected  for  the  future.  [Ref.l:p.3] 
The  goal  of  the  Navy  is  to  provide  for  standardized 
aircraft  and  engine  test  facilities  that  meet  both  safety  and 
environmental  concerns  while  providing  the  latest  technology 
to  monitor  the  engines.  To  accomplish  this,  complete 


development  of  test  facilities  with  reliable  operation  and 
maintenance  throughout  the  facility  life  cycle  is  necessary. 
[Ref .l:p.l] 


When  contemplating  the  inherent  design  criteria  for  the 
new  facilities,  problems  such  as:  turbulence,  non-concentric 
flow  patterns,  multiple  jet  engine  aircraft  (including  varying 
temperature  and  thrust  arrangements),  supersonic  velocities, 
extreme  temperature  and  pressure  gradients ,  and  awkward 
geometries  must  be  appraised.  [Ref. 5] 

Enigmas  associated  with  the  design  of  the  structure  may 
not  be  obvious.  Complete  information  dealing  with  the 
aerothermal  properties  is  a  must  in  basing  decisions  on  new 
multi-million  dollar  expenditures.  Therefore,  to  minimize  the 
possibility  of  an  erroneous  design,  the  PHOENICS  Code  is 
employed  first  on  documented  Test  Cell  and  Hush  House 
structures  having  actual  field  test  data.  After  verifying  the 
accuracy  of  PHOENICS  with  the  test  data,  the  Code  eventually 
could  be  used  at  a  Navy  installation  (such  as  the  Naval  Civil 
Engineering  Laboratory)  to  design  future  test  facilities 
providing  accurate  modeling  data. 

In  this  instance,  the  performance  of  a  Navy  J-79-GE-8  gas 
turbine  engine  located  at  Naval  Air  Station  Jacksonville, 
Florida  Hush  House  facility  is  investigated.  The  data  is  used 
to  determine  the  behavior  of  the  existing  structure  and  will 
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be  employed  to  assist  engineers  in  the  research  and 
development  (R&D)  phase  of  the  facility  design  project. 
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II .  PROBLiEayi  FORMULATION 


A.  HUSH  HOUSE  CHARACTERIZATION 

The  Hush  House  under  investigation  is  an  "Air  Force"  type 
located  at  Naval  Air  Station  Jacksonville,  Florida.  In  1980 
the  Air  Force  purchased  25  and  later  54  demountable  Hush 
Houses  with  Test  Cell  capability.  This  standardization 
accompanied  with  such  a  large  scale  purchase  and  the 
flexibility  of  aircraft  or  engine  test  capabilities  allowed 
for  a  cost  effective  combined  design.  Because  of  the  large 
scale  purchase  the  cost  per  Hush  House  was  a  relatively  low 
S2.2  million.  The  only  major  modification  performed  by  the 
Navy  to  the  Air  Force  design  was  inclusion  of  an  additional 
AFFF  fire  fighting  system.  The  Hush  House  design  does  not 
address  air  pollution  standards.  Some  type  of  air  pollution 
abatement  system  will  undoubtably  be  required  in  newer 
designs.  [Ref.l:p6] 

The  general  dimensions  of  the  Hush  House  is  28.7  meters 
(94  feet)  wide  by  25.5  meters  (83.7  feet)  long  with  the  oval 
shaped  "augmenter  tube"  extending  an  additional  31.3  meters 
(103  feet)  from  the  rear  of  the  Hush  House  (See  Figures  2.1 
and  2.2).  The  aircraft  is  placed  so  that  the  engine  exhaust 
is  approximately  4.6  meters  (15  feet)  from  the  augmenter  tube 
entrance.  The  most  common  construction  material  for  the 


structure  is  steel.  The  design  is  geometrically  simple 
allowing  for  mass  production.  The  air  flow  in  the  last  two 
inlets  nearest  the  augmenter  tube  provides  both  engine  intake 
and  cooling  air  for  the  surrounding  engine  area.  The  final  7 
meters  (23  feet)  of  augmenter  tube  is  angled  at  approximately 
45  degrees  running  from  the  floor  upward  to  provide  for  a 
smoother  directional  transition  of  gases.  Acoustic  damping 
materials  line  the  augmenter  tube,  inlets  and  aircraft  bay 
walls  to  reduce  exhaust  noise. 

The  Hush  House  is  used  in  support  of  the  Naval  Air  Rework 
Facility  (NARF)  located  at  the  Naval  Air  Station.  Although  the 
Hush  House  is  presently  utilized  for  engines  fully  mounted  in 
an  airframe,  it  is  anticipated  that  engine  "Test  Cell"  testing 
(where  an  engine  is  mounted  on  a  test  stand)  will  be  approved 
in  the  future.  The  facility  is  normally  in  operation  from  0700 
to  2400  Monday  through  Friday  with  an  average  of  1  to  5 
aircraft  being  tested  per  day.  One  of  the  most  common  aircraft 
tested  is  the  F/A-18  Hornet  powered  by  two  General  Electric 
F-404  gas  turbine  engines.  Typically,  after  completion  of 
overhaul  of  the  aircraft  and  engines,  the  aircraft  will  be 
towed  to  the  facility  and  backed  into  the  Hush  House.  The  hugh 
aircraft  inlet  doors  to  the  facility  will  then  be  closed  to 
provide  an  enclosed  area  free  of  foreign  matter,  outside  wind, 
rain,  and  other  environmental  and  personnel  concerns. 
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With  the  aircraft  fully  enclosed  and  securely  tied  down, 
a  normal  crew  of  three  mechanics  will  perform  final  engine 
"trim”.  One  person  is  located  in  the  Hush  House  Control  Booth 
to  ensure  manning  of  the  fire  fighting  systems.  The  second 
will  physically  adjust  the  fuel  control  mechanism  on  the 
engine  via  radio  signals  received  from  the  cockpit  where  the 
third  crewman  is.  The  engines  are  controlled  during  testing 
by  the  third  crewman  in  the  aircraft  who  is  additionally 
responsible  for  monitoring  the  various  fuel,  oil,  %  power, 
hydraulics,  and  other  gauges.  Testing  normally  lasts  for  a 
period  of  two  hours.  Only  one  engine  is  tested  at  a  time  with 
the  other  normally  at  idle.  Air  intake  and  cooling  is  provided 
by  five  inlet  areas  located  on  either  side  of  the  Hush  House 
(See  Figures  3.2  and  3.3).  Wire  screening,  air  diffuser  and 
noise  abatement  equipment  is  present  at  the  inlets.  Installed 
fire  fighting  and  wash-down  equipment  is  present  in  the 
aircraft  bay  area  with  the  air  inlets  capable  of  being 
automatically  closed. 

B.  ISSUES  CXDNCEaRNTNG  THE  HUSH  MOUSE 

The  Hush  House  is  currently  performing  as  designed  with 
the  below  exceptions.  Investigation  indicates  some  cracking 
of  the  welds  which  fuse  together  the  stainless  steel 
"corrugated"  sheet  metal  overlaying  the  sound  absorption 
material  located  in  the  augmenter  tube.  Also,  actuation  of 
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heat  sensors  in  the  overhead  of  the  aircraft  bay  sometimes 
occurs  when  the  engines  are  idled  for  a  period  of  time. 
Apparently  some  of  the  exhaust  gases  travel  upward  when  low 
mass  flow  rates  are  present  instead  of  going  down  the  31  plus 
meters  (103  feet)  of  augroenter  tube. 

In  general,  personnel  at  the  Naval  Air  Station  were 
pleased  with  the  overall  operation  of  the  Hush  House.  Although 
vibration  of  the  entire  physical  structure  has  been  noted  in 
this  facility  as  in  other  designs,  it  does  not  appear  to  cause 
any  structural  problems  such  as  the  requirement  for  inspection 
and  retightening  of  bolts  as  occurs  in  other  Navy  Hush  House 
designs . 

Physical  test  data  is  available  for  an  identical  Air  Force 
facility  for  the  F-4  Phantom  II  and  is  included  in  parts  of 
this  report  for  comparison  with  results  from  the  PHOENICS 
Code.  The  test  data  provides  an  important  comparison  for  this 
thesis  to  verify  the  computer  model.  [Ref. 2] 

A  physical  inspection  of  an  identical  Hush  House  to  Naval 
Air  Station  Jacksonville's  was  under-taken  in  June  1988  at 
Dobbins  Air  Force  Base  in  Atlanta  Georgia.  Discussions  with 
the  various  Air  Force  personnel  corroborated  the  previously 
mentioned  problems  with  the  exception  of  the  heat  sensor 
activation  problem.  The  Dobbins  Air  Force  Base  Hush  House  was 
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utilized  for  both  Test  Cell  and  Hush  House  operations 
involving  the  F-15  (EAGLE)  fighter. 
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Exhaust 


Figure  2.1  Actual 


III.  THE  PHCSaOXCS  CX>£>E 


A.  AN  ESCPLANATION 

PHOENICS,  which  stands  for  "Parabolic  Hyperbolic  Or 
Elliptic  Numerical  Integration  Code  Series",  is  a  computer 
code  capable  of  simulating  fluid  flow  (including  changes  of 
state),  heat  transfer  and  chemical  reactions.  The  code  is 
copyrighted  by  CHAM  (Concentration  Heat  and  Momentum) 
Limited  which  was  founded  by  D.  Brian  Spalding  in  1972. 
[Ref . 3 :p. 1 . 1] 

As  shown  in  Figure  3.1,  the  PHOENICS  Code  is  structured 
into  three  major  components  consisting  of  "EARTH" , 
"SATELLITE",  and  "GROUND  STATION".  The  "EARTH"  program  may  be 
considered  as  a  "blac)c  box"  which  contains  the  computational 
procedures  for  all  applications.  The  "SATELLITE"  program 
consists  of  several  data-input  methods  which  define  the 
boundaries  and  scope  of  the  problem.  Finally,  the  "GROUND 
STATION"  program  has  various  subroutines  containing  built- 
in  physical  models  which  may  or  may  not  be  selected  depending 
on  the  program  application. 

By  default  "EARTH"  can  solve  for  up  to  25  dependent 
variables.  In  the  Hush  House  application  these  include 
pressure,  all  three  dimensional  velocity  components, 
enthalpy  (temperature)  ,  turbulent  )cinetic  energy,  and  the 
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turbulent  kinetic  energy  dissipation  rate.  The  dependent 
variables  can  be  enlarged  if  desired. 


Figure  3.1  Satellite,  Earth  and  Ground 

The  PHOENICS  Code  is  exceptionally  flexible  allowing  for 
varying  sizes  and  scales  of  structures  from  "micro"  to 
"macro".  A  "final"  solution  is  highly  dependent  on  the 
accuracy  required  by  the  user  as  well  as  the  grid  sizing  used. 
The  higher  the  accuracy  and/or  number  of  "cells"  desired,  the 
greater  the  computational  effort  required  for  a  solution.  In 
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this  investigation,  a  "grid  independent  solution",  where  a 
sequence  of  "finer"  grids  would  be  used,  was  not  possible  due 
to  time  limitations. 

In  order  to  solve  a  given  problem  the  user  must  create  a 
"Ql”  file.  The  "Ql"  file  describes  all  aspects  of  the  system 
being  investigated.  This  file  is  partitioned  into  several 
groups  allowing  the  user  to  specify  various  aspects  of  the 
physical  system.  Examples  are  structure  size  and  shape. 
Additionally,  other  groups  are  capable  of  providing  the 
necessary  grid  spacing,  boundary  conditions,  engine  mass 
flow  rate,  engine  exhaust  temperature  and  printing  options. 
The  "Ql"  file  is  read  by  the  "SATELLITE"  program  which  is 
recognized  and  processed  by  "EARTH".  The  "EARTH"  program 
solves  discretized  conservation  equations.  A  typical 
conservation  equation  (before  discretization)  which  "EARTH" 
solves  can  be  expressed  in  the  form: 

5 

--  (rp<I>)  +  div(rpV$  -  rr»  grad  ♦)  =  rS» 

5t 

where : 

r  «  phase  volume  fraction 

♦  *  dependent  variable 

r*  ■  exchange  coefficient  (laminar  or  turbulent) 

S»  *  source  or  sink  term 
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p  *  density 
V  ■  velocity  vector 


Here  4  is  the  dependent  variable  currently  in  question.  For 
example,  using  the  conservation  of  energy  equation,  4  would 
be  enthalpy  and  r«  is  the  diffusion  constant.  [Ref . 3 :p. 2 . 4] 
The  values  computed  and  displayed  for  each  cell  are  the 
solutions  of  the  algebraic  equations  relevant  to  the  laws  of 
physics  for  that  cell.  The  laws  of  physics  are  those  of 
conservation  of  mass,  energy  and  momentum  and  any  other 
dependent  variable.  The  general  form  of  the  algebraic 
linearized  equations  resulting  from  a  discretization  of  the 
previous  field  equation  is: 


aE$i  +  aw^>y  +  Sn^n  +  as^s  +  an^B  +  ai^u  +  ar^T  +  S 


at  +  aw  +  aw  +  as  +  as  +  ai  +  ai  +  ap 


where  is  the  value  of  the  dependent  variable  at  a 

particular  cell  and  subscripts  E,  W,  N,  S,  H  and  L  are 
neighboring  locations.  The  values  of  ^  at  these  locations  are 
known  from  a  previous  iteration.  Quantities  S  and  ap  express 
the  influence  of  a  source  of  the  entity  4'.  The  a's  are 
influence  coefficients,  temporarily  treated  as  constants,  with 
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the  above  subscripts  expressing  interaction  between  neighbor 
cells.  [Ref . 3 :p. 2 . 4] 

Solving  the  multitude  of  algebraic  equations  is 
complicated  and  time  cons\iming.  An  iterative  manner  is  used 
until  the  imbalance  between  the  left  and  right  sides  of  all 
the  above  equations  are  negligible.  Several  schemes  can  be 
used  to  perform  the  "sweeps"  (where  a  sweep  solves  for  the 
dependent  variables  for  each  cell  in  the  grid  one  time) 
necessary  for  the  iteration  process.  The  results  will  be 
identical  although  the  time  to  complete  the  tas)t  will  vary 
with  the  iteration  scheme.  [Ref . 3 :p. 2 . 5] 

There  are  several  ways  to  determine  whether  the  problem 
has  "converged"  to  a  solution.  Results  are  monitored  to 
determine  whether  values  of  dependent  variables  are  remaining 
constant  or  fluctuating.  The  "error"  (residuals)  is  tabulated 
for  each  sweep  at  the  end  of  the  computer  printout.  By 
observing  "residuals"  from  each  sweep  a  pattern  of  decreasing 
("converging")  values  should  be  present.  If  the  "residuals" 
are  increasing  ("diverging")  adjustments  of  several  constants 
which  in-turn  apply  a  dampening  effect  upon  the  variables  in 
question  may  be  required.  The  printout  of  successive  residual 
values  demonstrate  convergence  by  a  reduction  in  the 
residuals.  That  is,  a  visually  downward  sloping  line  on  the 
accompanying  graph  of  residuals  shows  convergence. 
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Conservation  balances  on  each  dependent  variable  can  be 
verified  using  "net  source"  values  which  are  tabulated  at  the 
end  of  each  run.  The  mass  flow  rate  (mass  flow  rate  =  density 
*  area  *  velocity)  is  additionally  calculated,  in  this 
instance,  for  the  five  inlets  and  one  outlet  confirming  that 
no  hidden,  missing,  or  improper  boundaries  are  present.  The 
problem  has  "converged"  when  an  acceptable  level  of  error 
remains . 

Successful  "convergence"  of  the  PHOENICS  program  is 
dependent  upon  experience  and  an  understanding  of  the  results. 
If  the  residuals  of  a  dependent  variable  show  divergence,  the 
user  can  manipulate  the  rate  of  change  of  values  of  that 
variable  by  using  the  RELAX  command  to  return  the  variable 
back  tc  the  convergence  path.  Relaxation  is  placed  on  a 
variable  to  control  the  variables  maximum  change  from  each 
sweep.  This  minimizes  the  possibility  of  "over-correction" 
or  "blowing  up"  of  the  values.  Smaller  values  of  relaxation 
result  in  smaller  value  differences  per  iteration.  However, 
excessively  small  relaxation  values  require  additional 
computer  time.  A  balance  gained  by  user  experience  enables 
optimal  convergence  with  minimal  computer  time. 

When  multiple  runs  are  required  for  large  programs  output 
data  is  stored  in  a  "DF09"  file  which  allows  the  next  run  to 
resume  the  computations  (i.e.,  a  restart)  where  they  last 
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ended.  The  complexity  of  the  problem  and  the  required 
accuracy  determine  the  quantity  of  computer  runs  needed.  The 
number  of  runs  can  vary  from  the  tens  to  the  thousands  and 
higher . 

B-  ASSUMPTIONS  IN  THE  MOBEIj 

There  are  several  approximations  that  deserve  mention  in 
the  modeling  of  the  Naval  Air  Station,  Jacksonville,  Florida 
Hush  House.  Because  the  geometry  of  the  model  uses  cartesian 
coordinates,  the  actual  Hush  House's  rounded  and  curved 
surfaces,  including  the  augmenter  tube.  Hush  House  overhead, 
engine  and  exhaust  gas  exit  ramp,  were  approximated  as  flat 
stepped  surfaces. 

The  inlet  "muffler"  arrangement  has  been  eliminated  for 
simplification.  The  corrugated  metal  found  throughout  the 
structure  to  cover  the  sound  absorption  material  was 
considered  to  have  a  friction  surface.  The  model,  therefore, 
included  friction  for  all  walls,  floors,  etc. 

The  actual  test  data  used  for  comparison  between  the 
model  and  facility  contain  several  notable  differences.  The 
model  uses  the  symmetry  of  the  centerline  of  the  Hush  House 
to  reduce  the  number  of  cells  by  half.  The  engine  is  located 
on  the  centerline  in  an  ideal  cradle  position.  Because  of 
the  symmetry  only  half  of  the  engine  is  observed  as  if  the 
Hush  House  was  being  used  as  a  test  cell  without  the  aircraft 


18 


present.  Since  the  correct  engine  exhaust  area  is  used  the 
units  of  mass  flow  rate  and  temperature  are  correct  for  the 
model  configuration.  The  engine  is  configured  for  afterburner 
operation.  The  physical  aircraft  is  not  present. 

Actual  test  data  used  for  comparison  consists  of  the 
entire  aircraft  being  placed  in  the  facility.  The  engine 
configuration  has  both  engines  running  with  one  at 
afterburner  while  the  other  engine  is  at  "engine  idle".  The 
obvious  result  of  the  engine  arrangement  is  a  non-symmetrical 
engine  location  in  comparison  to  the  symmetric  model. 

C.  MODELiXNG  THE  PRC^LEM 

Data,  including  blueprints  of  the  Hush  House  and  aircraft 
engine  information,  was  gathered  to  accurately  depict  the 
facility  configuration.  The  information  obtained  was  in 
English  units  and  was  converted  to  the  International  Metric 
System  or  "S.I."  using  standard  conversion  factors. 

With  the  appropriate  information  gathered  it  was  time  to 
"fit"  the  Hush  House  to  the  cartesian  coordinate  system  of  X, 
Y  and  Z  coordinates  using  PHOENICS  version  1.3.  The  cartesian 
coordinate  system  was  choosen  because  of  its  simplicity.  Since 
the  actual  structure  does  have  rounded  surfaces  some 
approximations  were  made  specifically  to  the  augmenter  tube 
and  overhead  of  the  inlet  bay.  Additionally,  the  aircraft  was 
eliminated  to  reduce  the  complexity  of  the  problem  allowing 
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the  structure  to  serve  in  its  "Test  Cell"  configuration. 
Drawings  are  found  in  Figures  3.2  and  3.3  of  this  section 
illustrating  the  key  items  found  in  the  structure.  The 
drawings  view  through  the  structure  to  note  the  key  items, 
such  as:  air  inlets,  augmenter  tube  and  engine  locations. 
Additionally,  the  structure  walls  were  considered  to  be  of 
infinitesimal  thickness  to  further  simplify  the  grid 
structure.  The  engine  was  simplified  from  a  cylindrical  shape 
to  one  of  a  rectangle  of  equal  intake  and  exhaust  area  as  the 
original  cylinder.  Other  small  items,  such  as  the  engine 
stand,  were  eliminated. 

To  take  advantage  of  the  symmetry  of  the  Hush  House  and 
thereby  reduce  computer  CPU  time,  only  half  of  the  Hush  House, 
engine  and  augmenter  tube  were  modelled  in  the  "Ql"  file.  This 
was  easily  accomplished  by  dividing  the  structure  at  the 
centerline.  This  also  appreciably  reduces  the  size  of  the 
"Ql"  file. 

Figures  3.2  and  3.3  graphically  display  the  X-Z  and  the 
y-Z  planes  of  the  structure.  Intuition,  trial-and-error ,  and 
available  computer  time  dictated  the  grid  spacing  used. 
Obviously,  volume  areas  of  interest  such  as  the  engine  exhaust 
and  augmenter  required  smaller  cell  size. 

The  "Ql"  file  describes  the  structure  using  several  types 
of  statements.  The  "CONPOR"  statement  allows  the  user  to 
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define  walls,  overheads  and  the  floor.  The  "PATCH"  and 
"COVAL"  statements  work  together  to  indicate  items  such  as 
an  engine  exhaust  or  atmospheric  boundaries.  Dependent 
variable  values  for  say  the  atmospheric  boundaries  can  be 
"fixed"  or  allowed  to  "float"  depending  on  the  problem 
definition.  When  the  walls,  inlet  bay  overhead,  augmenter 
tube,  engine  and  floor  dimensions  have  been  measured,  a  grid 
for  the  system  can  be  constructed,  based  on  locations  and 
anticipated  responses  over  different  "areas"  of  the  system. 
Areas  of  anticipated  large  response  gradients  would  have 
smaller  cell  volumes  than  low  activity  areas. 

Note  that  in  figure  3.2  the  Z-axis  consists  of  40  cells 
consecutively  numbered  from  1  to  40  going  left  to  right.  The 
X-axis  consists  of  10  cells  consecutively  numbered  1  through 
10  going  from  lower  to  upper.  Figure  3.3  is  similarly  numbered 
where  the  Y-axis  runs  from  1  to  12  going  from  lower  to  upper 
and  the  Z-a:-;is  is  the  same  as  previously  described. 

In  addition  to  describing  the  geometric  structure  and 
prescriting  the  "P.ELAXATION"  devices,  the  properties  of  the 
fluid  medium  must  be  specified.  The  following  boundary 
conditions  were  used; 

1.  ambient  temperature  =  295  K 

2.  am.bient  pressure  =  101325  Pa 

3.  gas  constant  =  286  Joules/ (kg  -  K) 


4.  Cp  *  1004  Watt-sec/ (kg  -  K) 

5.  mass  flow  rate  (engine)  =  197.72  kg/(sec-sq  meters) 

6.  temperature  (engine  out)  »  1944  K 

Additionally,  "friction**  was  applied  to  all  surfaces  of 
the  Hush  House.  Appendix  A  contains  the  final  form  of  the 
Q1  file  for  reference.  An  in-depth  explanation  of  all  the  "Ql" 
program  parameters  can  be  found  in  the  "Beginners  Guide  and 
Reference  Manual"  (See  Reference  4)  . 

Several  methods  and  theories  have  been  previously  used  to 
attempt  to  reach  a  "converged"  solution  in  other  test 
facilities.  The  method  chosen  for  this  project  was  to 
initially  minimize  the  "excitation"  of  the  model  thereby 
reducing  the  "jump"  required  in  solving  the  momentum, 
enthalpy,  turbulent  kinetic  energy  dissipation  rate,  turbulent 
kinetic  energy  and  mass  equations.  To  accomplish  this,  during 
successive  runs  (where  a  "run"  consists  of  several  "sweeps"), 
small  incremental  increases  to  the  source  were  applied,  i .  e., 
jet  engine  exhaust  temperature  and  mass  flow  rate  were  slowly 
increased,  from  run  to  run,  while  carefully  monitoring  the 
residuals  as  well  as  the  mass  flow  rate  in  and  out  of  the 
structure.  The  turbulence  model,  "KEMODL" ,  was  used 
throughout  rather  than  considering  the  flow  first  as  laminar 
and  then  adding  turbulence.  This  method  appears  to  have  met 
with  considerable  success.  The  table  that  follows  illustrates 
the  iterative  process  used  throughout  the  computer  modeling: 
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Table  3.1  ITERATIVE  PROCESS 


Temperature 

1  Flow  Rate 

1  Sweeps 

400 

i 

j  600 

500 

1 

1  600 

750 

1 

1  600 

750 

1  50 

1  600 

1000 

1 

1  2500 

1250 

1  75 

1  1200 

1250 

1  100 

1  1200 

1250 

1  125 

1  1200 

1500 

1  125 

1  1200 

1750 

1  125 

1  2000 

1750 

1  150 

1  2000 

1944 

J  150 

1  2000 

1944 

1  l‘^5 

1  5000 

1944 

1  197.72 

1  5500 

^  I  I.  I.  ■  ■  - 

Total  Sweeps  =  26,200 


This  amounts  to  appro-ximately  147  hours  of  computer  CPU  time 
on  the  IBM  3033  mainframe  computer  installed  at  NPS. 

Using  the  MVS  "batch"  system  on  the  IBM  3033  mainframe 
proved  to  be  time  consuming.  The  PHOENICS  Code  is  designed  to 
work  both  interactively  and  in  batch  mode.  Due  to  the  high 
computer  usage  at  NPS  the  "interactive"  mode  was  not 
available.  Limited  to  the  batch  mode,  careful  planning  for 
evening  or  night  runs  was  required  to  maximize  CPU  time.  The 
number  of  "Degrees  of  Freedom"  (DOF)  that  must  be  solved  for 
gives  an  idea  of  the  problem  complexity.  DOF  =  NX  *  NY  *  NZ 
*  the  7  dependent  variables.  Solving  this  for  the  grid  size 
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of  10  *  12  *  41  yields  nearly  35,000  values  calculated  per 
each  sweep.  Typically,  150  to  200  sweeps  could  be  accomplished 
during  the  day  using  four  15  minute  (class  "G")  runs  of  50 
sweeps  each.  After  normal  working  hours  three  one  hour  (class 
"J")  runs  totalling  500  to  600  sweeps  were  accomplished.  Data 
stored  on  the  all  important  DF09  file  was  continuously 
"backed-up”  after  a  successful  run  to  prevent  loss  of  data  if 
an  improper  subsequent  run  resulted  in  erroneous  data.  The 
backed-up  data  could  then  be  restored  to  the  DF09  file  to 
enable  continuation  of  the  runs  without  starting  at  the  very 
beginning. 
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Figure  3.2  Hush  House  Model  (X-Z  plane) 
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IV.  REISULTS 


A.  THE  FINAL  SOXJUTIOM 

After  obtaining  the  "solution”  for  the  problem  the 
question  still  remaining  is  just  how  "good”  are  the  computer 
model  results?  Some  of  the  "tools"  available  for  use,  such  as 
checks  on  residual  values  and  mass  flow  rate  calculations, 
were  discussed  in  Chapter  III.  The  checks  indicate  that  a 
converged  solution  is  obtained,  but,  still  do  not  prove  the 
computer  model  data  validity  with  the  actual  facility  data. 

Using  actual  (experimental  field  measurements)  results 
from  Reference  2  allows  for  an  accurate  evaluation  of  the 
computer  model.  Appendices  ”C"  and  "D"  offer  some  experimental 
data  for  velocities  and  temperatures,  thereby  indicating  mass 
flow  rate  and  augmentation  ratio.  Additionally,  figures  4.1, 
4.2  and  4,3  found  in  this  section  give  the  numerical 
(computer)  model  results  of  differential  pressure,  temperature 
and  axial  velocity  profiles  encountered  at  the  IX=1  slab. 
Graphics  generated  pictures  (generated  by  "PHOTON"  which 
converts  the  computer  generated  results  to  graphical  format) 
allow  for  the  interpretation  and  digestion  of  the  nearly  300 
pages  of  numerical  output  of  a  PHOENICS  run.  A  condensed 
output  is  found  in  Appendix  "B".  In  this  case,  PHOTON  graphics 
specifically  provides  graphical  information  for  a  better 
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interpretation  of  velocity  fields  and  both  pressure  and 
temperature  profiles  and  contours  enabling  a  better 
understanding  of  computer  generated  data.  Use  of  all  these 
additional  "tools"  allow  for  a  determination  as  to  whether  the 
computer  results  are  valid. 

Observing  figures  4.4  and  4.5  (top  looking  down)  show  the 
flow  pattern  of  the  air  entering  the  five  air  inlets  into  the 
Hush  House  at  IY=4  (engine  level) .  k  "typical"  low  velocity 
air  flow  region  (Figure  4.4),  illustrated  by  the  minuscule 
length  of  the  velocity  arrows  (velocities  of  less  than  2 
meters/sec  or  6.5  feet/sec),  is  noted  between  the  first  two 
and  final  three  inlets.  Figures  4.6,  4.7  and  4.8  illustrate 
the  side  view  velocity  flow  patterns  for  the  overall  facility, 
engine  area,  and  exhaust  to  atmosphere  respectively.  Note  the 
ever  reducing  velocity  magnitudes  which  vary  from  818  meters 
per  second  (2684  feet/sec)  at  the  engine  exhaust  to  less  than 
100  meters  per  second  (328  feet/sec)  as  the  gases  travel 
further  away  from  the  engine  and  down  the  augmenter  tube  until 
reaching  the  exhaust. 

Figures  4.4,  4.5  and  4.9  (top  looking  down)  show  the  air 
and  gas  flow  patterns  entering  and  exiting  the  engine  region 
and  augmenter  tube  at  engine  level.  Figure  4.10  indicates  the 
flow  patterns  at  the  top  of  the  augmenter  tube  (iy=8)  for  the 
entire  Hush  House.  Of  interest  is  the  counter-flowing  currents 
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at  the  top  of  the  augmenter  tube  starting  at  the  engine 
exhaust  and  running  to  nearly  half  way  down  the  augmenter  tube 
itself . 

The  engine  occupies  cells  defined  by  IX«1,  IY«4  to  5  and 
IZsl2  to  14.  As  anticipated,  figures  4.4,  4.9  and  4.11  show 
that  at  engine  level  the  air  in  the  Hush  House  moves  toward 
the  engine  inlet  from  the  surrounding  area  while  air  also  is 
entrained  at  the  exhaust  of  the  engine.  Figure  4.1  illustrates 
this  entrained  "effect"  noted  by  the  high  differential 
pressure  of  the  exhaust  gases  immediately  exiting  the  engine 
at  IX=1  and  IY=4.  Less  than  a  meter  after  exiting  the  engine 
exhaust,  the  exhaust  gases  and  inlet  air  have  begun  mixing 
which  results  in  an  area  of  drastic  temperature  reduction  for 
the  exhaust  gases  in  the  Hush  House  and  immediate  augmenter 
area.  Gradually,  as  the  gases  pass  further  downstream, 
pressure  begins  rising  which  reduces  the  velocity  of  the  gases 
in  the  augmenter  tube  until  the  gases  exit  to  the  atmosphere. 
At  the  inlet  area  to  the  engine  (at  Z=18  meters  on  figure  4.1) 
a  low  pressure  area  is  noted  as  one  would  also  anticipate  as 
the  air  is  drawn  into  the  engine  inlet. 

The  exhaust  gases  exit  the  engine  at  1944  Kelvin  (3040  F) . 
The  inlet  air  and  exhaust  gases  mix  due  to  the  entrainment 
"effect”  causing  a  distinct  drop  (See  figure  4.2)  in 
temperature  as  the  gases  travel  further  down  the  augmenter 
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tube.  Temperature  contours  at  the  IX=1  slab  shown  in  figures 


4.12  and  4.13  portray  the  temperature  patterns  out  of  the 
engine  progressing  down  and  out  of  the  augmenter  tube. 
Analogous  temperature  contours  at  the  IXb3  slab  are  presented 
in  figures  4.14  and  4.15.  Figures  4.16,  4.17  and  4.18  embody 
the  temperature  contours  at  the  lY^s  (engine)  level  while 
figure  4.19  shows  the  temperatures  at  the  top  of  the  augmenter 
tube  at  IY=8.  The  temperature  contours  indicate  a 
progressively  widening  area  of  decreasing  exhaust  gas 
temperatures  where  the  vast  majority  of  the  temperature  drop 
occurs  in  the  first  five  meters  after  the  engine  exhaust  (the 
temperature  differential  goes  from  1944  K  to  800  K) . 

Differential  pressure  contours  of  the  engine  exhaust  at 
IX=1  are  shown  in  figures  4.20  (IZ=  14  to  24)  and  4.21  (IZ® 
25  to  40).  As  indicated  in  the  pressure  profile,  the  exhaust 
jumps  to  nearly  6000  pascals  above  ambient  upon  exiting  the 
engine,  but  soon  after  drops  into  a  low  pressure  region  until 
reaching  ambient  pressure  half  way  down  the  augmenter  tube  at 
IZ=24,  Gradually,  the  pressure  builds  until  the  gases  exhaust 
to  the  atmosphere.  Similar  contours  for  the  IX=3  slab  are 
included  in  figures  4.22  and  4.23.  Pressure  contours  for  the 
X-Z  planes  at  IY=5  (engine  level)  and  IY=8  (top  of  the 
augmenter  tube)  are  contained  in  figures  4.24,  4.25  and  4.26. 
Of  noted  interest  is  how  the  differential  pressures  remain 
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fairly  uniform  at  a  given  IZ  location  as  the  gases  travel  down 
the  augmenter  tube  with  the  exception  of  the  region  initially 
downstream  from  the  engine  exhaust. 

B.  ACTUAL  VERSUS  M3DEL.  OCMPARZSCTTS 

Generally,  the  results  from  the  computer  model  closely 
approximated  the  actual.  The  below  "actual"  data,  and  that  of 
Appendix  "C"  and  "D"  for  the  computer  model  are  indicated  in 
both  "English"  and  "S.I."  units  below.  A  comparison  of  the 
PHOENICS  output  visually  shown  in  Appendices  "C"  and  "D" 
illustrates  the  dissimilarity  between  the  actual  and  computer 
model.  As  indicated  by  the  data,  the  differences  in  the 
computer  model  (235.5  ft/sec  or  71.8  meters/sec)  and  actual 
(227.2  ft/sec  or  69.3  meters/sec)  averaged  velocities  found 
in  the  augmenter  tube  at  slab  IZ=27  differed  by  less  than  4 
percent.  Averaged  temperatures  in  the  same  location  differed 
between  401.6  degrees  Fahrenheit  (478  K)  for  the  computer 
model  and  454  degrees  Fahrenheit  (507  K)  actual  with  a 
resultant  difference  of  6  percent.  Reference  2  does  not 
indicate  whether  the  effects  of  radiant  energy  from  the  engine 
exhaust  are  considered  for  the  actual  temperature  readings 
taken  in  the  augmenter  tube.  If  these  effects  were  not 
considered  the  6  percent  difference  would  be  reduced.  The  Hush 
House  mass  flow  rate  for  the  computer  model  (2373  Ibm/sec  or 
1076  kg/sec)  versus  the  actual  (2067  Ibm/sec  or  938  kg/sec) 


31 


differed  by  14.8  percent.  The  computer  model  mass  flow  rate 
value  was  actually  doubled  (hence  the  2373  Ibm/sec  or  1076 
kg/sec  value)  to  "model"  the  entire  engine  and  Hush  House. 
The  Augmentation  Ratio  is  defined  as  follows: 


A.R. 


M 1 D 1 e  t  ~  Me  n  g i n  e 


Ml  Diet 


where : 

Miniet  =  inlet  mass  flow  rate 

Mengine  =  engine  mass  flow  rate 

The  inlet  mass  flow  rate  =  air  density  *  five  inlet 
opening  areas  *  averaged  inlet  velocity.  The  engine  mass  flow 
rate  is  given  at  197.72  kg/sec/sq  meters  *  0.195  sq  meters 
engine  exhaust  area  for  half  of  the  engine. 

Augmentation  ratios  range  from  5  to  15  depending  on  the 
design  of  the  Hush  House  or  Test  Cell  where  the  higher  the 
ratio  typically  indicates  a  more  efficient  design.  The 
computer  model  augmentation  ratio  was  found  to  be  10.6  while 
the  actual  was  8.3.  The  resultant  difference  is  therefore 
nearly  28  percent.  Considering  the  limited  data  taken  at  this 
actual  power  level  (See  Appendix  C)  it  is  quite  conceivable 
that  the  actual  augmentation  ratio  could  be  greater. 


For  the  actual  results  the  velocity,  pressure  and 
temperature  measuring  devices  were  photographed  because  of  the 


large  fluctuations  of  the  gauges.  Additionally,  velocity  data 
at  the  final  three  inlets  could  not  accurately  be  taken  making 
it  difficult  to  accurately  determine  the  augmentation  ratio 
and  the  true  velocities  at  the  inlets. 


Figure  4.1  Differential  Pressure  Profile  at  Volume  Slab 
IX^l  for  1Y«4  (Engine  Height) 
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NPS  JACKSONVILLE  HUSH  CHAMBER 


Figure  4.4  Flow  Vector  Plot  of  Hush  House  Inlets  and 
Engine  Area  {X-Z  Plane  at  IY=4) 
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Figure  4.5  Flow  Vector  Plot  in  X-Z  Plane  at  IYs4 
(Engine  Exhaust  and  Augmenter  Tube) 
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NftS  JPCKSONVILLE  HUSH  CHAMBER  PHOENICS 


Figure  4.6  Flow  Vector  Plot  in  Y-Z  Plane  at  IX**! 

(Entire  Grid) 


39 


NftS  JACKSONVILLE  HUSH  CHAMBER  PHOENICS 


VELOCITY  VECTORS  IN  Y-Z  PLPNE  PT  IX 


N 


Figure  4.7  Flow  Vector  Plot  in  Y-Z  Plane  at  IX*1 
(Engine  Exhaust  Area) 
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Figure  4.8  Flow  Vector  Plot  in  Y-Z  Plane  at  IX>1 
(Exhaust  to  Atmosphere) 
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Figure  4.10  Flow  Vector  Plot  in  X-Z  Plane  at  iy=8 
(Mirrored  Image  of  Entire  Grid) 
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Figure  4.11  Flow  Vector  Plot  of  Augmenter  Tube  Area 
for  X-Z  Plane  at  iy-4 
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Figure  4.12  Temperature  Contour  in  Y-Z  Plane  at  IX«1 
(Engine  and  Augment er  Tube) 


NAS  JACKSONVILLE  HUSH  CHAMBER  PHOENICS 


Figure  4.13  Temperature  Contour  in  Y-Z  Plane  at  IX«=1 
(Augmenter  Tube  and  Exit) 
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NAS  JACKSONVILLE  HUSH  CHAMBER  PHOENICS 


Figure  4.14  Temperature  Contour  in  Y-Z  Plane  at  IX*3 

(Engine  Exhaust) 
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NPS  JACKSONVILLE  HUSH  CHAMBER  PHOENICS 


Figure  4.15  Temperature  Contour  in  Y-Z  Plane  at  IX=3 
(Engine  Exhaust  to  Atmosphere) 
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NftS  jpcksonville  hush  chamber  phoenics 


Figure  4.17  Temperature  Contour  in  X-Z  Plane  at  IY»5 

(Engine  Exhaust) 
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TEMPERPTURE  CONTORS  IN  X-Z  PLPNE  PT  lY 


Figure  4.18  Temperature  Contour  in  X-Z  Plane  at  IY*5 
(Engine  Exhaust  in  Augmenter  Tube) 
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NPS  JPCKSONVILLE  HUSH  CHPfflBER  PHOENICS 


Figure  4.19  Temperature  Contour  in  X-Z  Plane  at  IY=8 
(Engine  Exhaust  in  Augmenter  Tube) 


52 


NfiS  JACKSONVILLE  HUSH  CHAMBER  PHOENICS 


Figure  4.20  Differential  Pressure  Contour  in  Y-Z  Plane 

at  IX=1  (Engine  Exhaust) 
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NftS  JPCKSONVILLE  HUSH  CHPHBER  PHOENICS 


Figure  4.22  Differential  Pressure  Contour  in  Y-Z  Plane 
at  IX=3  (Engine  Exhaust  and  Augmenter  Tube) 
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NPS  JACKSONVILLE  HUSH  CHO^SER  PHOENICS 


Figure  4.23  Differential  Pressure  Contour  in  Y-Z  Plane 
at  IX=3  (Exhaust  to  Atmosphere) 
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NAS  JACKSONVILLE  HUSH  CHAMBER  PHOENICS 


Figure  4.25  Differential  Pressure  Contour  in  X-Z  Plane 
at  iy=5  (Auginenter  Tube) 
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^IftS  JPCK50NVILLE  HUSH  CHOM8ER  PHOENICS 


Figure  4.26  Differential  Pressure  Contour  in  X-Z  Plane 
at  IY=8  {Augmenter  Tube) 
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NOS  JACKSONVILLE  HUSH  CHAMBER  PHOENICS 


V.  SUMlyD\RV  AND  REXXMMEasiDATXONS 


A.  SUMMARY 

Using  the  Parabolic,  Hyperbolic,  or  Elliptic  Numerical 
Integration  Code  Series  (PHOENICS)  a  computer  model  of  the 
Naval  Air  Station  Jacksonville  Florida's  Hush  House  facility 
with  a  J-79-GE-8  gas  turbine  engine  was  simulated.  The 
performance  of  this  computer  model  was  compared  to  actual 
performance  data  taken  with  an  F-4  Phantom  II  undergoing 
testing  at  an  identical  facility  located  at  Kelly  AFB.  An  on¬ 
site  inspection  was  conducted  at  Dobbins  AFB  in  Atlanta 
Georgia  to  provide  "hands-on"  experience  with  the  facility. 

The  computer  analysis  consisted  of  slowly  building  up 
excitation  in  the  model  (using  the  KEMODL  turbulence  model) 
until  full  power  was  obtained.  Computer  generated  (tabulated 
and  graphical)  results  were  compared  with  the  actual  test 
data.  Results  of  the  computer  model,  in  comparison  to  the 
actual  data,  were  in  good  agreement.  Flow  directions, 
velocities  and  temperatures  were  consistent  with  actual 
results . 

B  -  R£XX3MME3«]DA'rxahlS 

A  comparison  of  actual  to  computer  model  data  for  this 
thesis  indicates  that  PHOENICS  is  a  useful  design  and 


troubleshooting  tool  allowing  the  user  to  better  understand 
the  physics  of  the  model.  PHOEMICS  at  NFS  is  now  just  coming 
into  potent  use  with  a  variety  of  research  projects  for 
aeronautical  and  mechanical  engineers. 

To  ensure  proper  and  full  use  of  this  code  several  major 
"glitches"  must  be  overcome.  Time  constraints  placed  on  the 
use  of  the  IBM  3033  main-frame  computer  limit  the  number  of 
"sweeps"  obtainable  by  the  student.  The  batch  system  does  not 
lend  itself  to  correcting  possible  errors  that  may  arise.  An 
interactive  system,  such  as  the  "VAX"  stand  alone  systems 
being  installed  in  Mechanical  Engineering,  is  an  ideal  answer. 

PHOTON  was  not  available  with  the  PHOENICS  installation 
on  the  IBM  3033  main-frame  computer.  A  complicated  conversion 
process  (See  Appendix  "E")  was  undertaken  to  convert  IBM  to 
VAX  "format"  to  allow  the  use  of  PHOTON  using  data  from  the 
main-frame.  This  conversion  process  which  involved  a  swap  of 
information  from  IBM  to  DEC  to  VAX  hardware  took  nearly  2 
weeks  to  accomplish.  In  September  1988  this  transfer  procedure 
will  no  longer  be  possible  at  NPS  with  the  elimination  of 
certain  hardware  at  the  Computer  Science  Department  although 
data  transfer  will  still  be  possible  to  other  facilities  such 
as  CHAM. 

The  use  of  PHOTON  early  in  the  computer  "runs"  would  have 
uncovered  several  missing  "CONPOR"  statements  and  therefore 
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"verified"  the  structural  integrity  of  the  facility. 
Additionally,  the  amount  of  time  required  to  "wade"  through 
the  300  pages  of  printout  could  be  reduced  to  several  video 
screens  with  a  much  reduced  volume  of  printout. 

The  computer  model  used  in  this  thesis  was  formulated 
using  cartesian  coordinates.  Several  parts  of  the  structure, 
including  the  augmenter  tube  and  engine,  were  curve  shaped. 
Geometric  approximations  were  made  to  substitute  for  the  curve 
shaped  areas.  A  more  computationally  time  consuming  model  with 
increased  numbers  of  cells  and/or  the  use  of  body-fitted- 
coordinates  (BFC)  would  enhance  the  computer  model. 

As  noted  in  Figure  4.2,  over  fifty  percent  of  the 
reduction  in  engine  exhaust  gas  temperatures  occurred  in  the 
first  20  percent  of  the  augmenter  tube.  Could  an  improved 
design  include  shortening  the  augmenter  tube?  This  type  of 
"what  if"  question  can  be  answered  using  PHOENICS.  With  the 
variable  data  file  ("DF09")  already  in  place  the  grid  pattern 
can  be  modified  to  shorten  the  augmenter  tube.  Minimal 
computations  would  be  required  to  gather  definitive  results. 

Parameter  studies,  including  shortening  the  length  of  the 
augmenter  tube  (to  determine  the  effect  of  augmenter  tube 
length  on  performance) ,  movement  of  the  engine  location  in  the 
Hush  House,  changing  the  inlet  locations  and  sizes,  and  the 
use  of  different  turbulence  models  could  lead  to  improvements 
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on  the  present  Hush  House  design.  In  fact,  the  ability  to 
change  parameters  and  observe  their  effect  on  system 
performance  is  the  prime  attribute  of  numerical  modeling. 

Finally,  a  discussion  of  turbulence  is  in  order.  PHOENICS 
offers  four  turbulence  models  which  can  be  readily  selected. 
The  four  models  are:  Algebraic,  Mixing  Length,  K-Mixing  Length 
or  K-Epsilon.  The  user  can  also  formulate  a  model  of  his  own 
choosing  which  can  be  incorporated  into  the  code  via  the 
"GROUND"  program.  Two  questions  arise.  First,  which  turbulence 
model  is  appropriate  for  the  problem  under  investigation? 
Secondly,  how  strong  is  the  coupling  between  other  system 
response  variables  and  turbulence  variables?  The  second 
question  might  be  answered  by  running  the  PHOENICS  Code  on  the 
problem  using  several  different  turbulence  models. 
Unfortunately  there  was  not  adequate  time  to  perform  this  in 
the  investigation.  The  first  question  can  not  be  resolved 
until  a  better  understanding  of  turbulence  is  obtained. 

C-  CX3NC1AJSIC3NS 

The  computer  model  results  generated  by  this  thesis 
closely  parallel  the  actual  test  data.  The  results  therefore 
indicate  that  the  PHOENICS  code  can  be  successfully  utilized 
to  model  future  Test  Cell  and  Hush  House  applications  as  well 
as  "trouble-shoot"  suspect  areas  of  interest.  Also,  areas  that 
may  be  difficult  to  properly  monitor  with  test  equipment  could 


be  modelled  instead.  With  the  use  of  additional  "cells"  and/or 
BFC  a  more  refined  model  could  more  closely  duplicate  the  test 
facility  providing  more  accurate  results.  Follow-up  parametric 
studies  to  improve  and  verify  the  design  and  data  seem  in 
order.  The  results  indicate  that  the  Hush  House  is  performing 
as  intended. 
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APPEIODXX  A 
THE  PHOENICS  Ql  FTLE 


TALXsF,RUN(l,l) 


»  NAVAL  AIR  STATION  JACKSONVILLE  FLORIDA  HUSH  HOUSE  TEST 

•  FACILITY.  THESIS  PROJECT  OF  LT  ERIC  A.  NICOLAUS  USCG 

•  HITH  THESIS  INSTRUCTOR  PROF  DAVID  SALINAS. 


«  LIST  OF  VARIABLES 

» 

•  HJET:  VELOCITY  OF  ENGINE  EXHAUST  IMETERS/SEC) 

»  MJET:  MASS  FLOH  RATE  OF  GAS  I KGw/SECOND/SQ  METER) 

•  TJET:  JET  EXHAUST  GAS  TEMPERATURE  I  DEGREES  KELVIN) 

•  PJET:  EXHAUST  PRESSURE  OF  JET  ENGINE  INEHTONS/SQ  METER) 

•  RHOJET:  DENSITY  OF  THE  GAS  (KG/CUBIC  METER) 

•  HJET:  ENTHALPY  I  JOULES/KG  OR  NEHTON«nETERS/KC ) 

«  KEJET:  KINETIC  ENERGY  (JOULES  OR  NEI4TON-METERS ) 

•  GASCON:  GAS  CONSTANT  ( JOULES/I KG-DEG  KELVIN)) 

«  TAMB:  AMBIENT  AIR  TEMPERATURE  (DECREES  KELVIN) 

•  HAKB:  AMBIENT  ENTHALPY  (J0ULES/X6) 

•  KEINJ:  TURBULENT  KINETIC  ENERGY  (JOULES)  OF  JET 

«  CSUBP;  SP  HT  S  CONST  PRE  ( JOULES/( KG-DEG  KELVIN)) 

•  EPINJ:  DISSIPATION  RATE  OF  TURBULENT  KE  (JOULES)  OF  JET 
«  KEINA:  TURBULENT  KINETIC  ENERGY  (JOULES)  AMBIENT 

«  EPINA:  DISSIPATION  RATE  OF  TURBULENT  KE  (JOULES)  AMB 
«  PI:  PRESSURE  ( NEWTONS/SOUARE  METER) 

*  HI:  ENTHALPY  < JOULES/KG) 

*  Ul:  VELOCITY  IN  X-DIRECHON  (METERS/SEC) 

*  VI:  VELOCITY  IN  Y-OIRECTION  (METERS/SEC) 

*  HI:  VELOCITY  IN  Z-DIRECHON  (METERS/SEC) 

«  TMPl:  TEMP  OF  THE  FIRST  PHASE  (DEGREES  KELVIN) 

*  RHOl:  DENSITY  OF  THE  FIRST  PHASE  (KG/CUBIC  METER  I 

*  PRESSO:  REFERENCE  ATMOSPHERIC  PRESSURE  (NEWTONS/SG  mir) 
«  TMPIA:  GROUND  REF  TEMP  USED  IN  GRNDZ  (DEG  KELVIN) 

»  TMPIB:  1/CSUBP  USED  IN  GRF(0Z  EQU  IKC-DE6  KELVIN/JOULES) 
«  DRHIOP:  DELTA  RHO/OELTA  PRESSURE  H/RESPECT  TO  HI 

*  FACT:  CORRECTION  OR  "FUDGE"  FACTOR 

*  ENUL:  LAMINAR  KINEMATIC  VISCOSITY  (METERS  SQ/SEC ) 

*  ENUT:  TURBULENT  KINEMATIC  VISCOSITY  (METERS  SQ/SEC) 


GROUP  I.  RUN  TITLE  AND  OTHER  PRELIMINARIES 
TEXKNAS  JACKSONVILLE  HUSH  CHAMBER) 

REAL(  HJET, MJET, TJET. PJET, RHOJET,HJET,KEJET>GASCON>CSUBP, TAMB) 
REAL(HAM&,KEINJ,EPINJ,KEINA,EPINA) 

GROUP  2.  TRANSIENCE)  TIME-STEP  SPECIFICATION  ( OEFAULT*STEADY ) 

GROUP  3.  X-DIRECTION  GRID  SPECIFICATION 
NX^IO 

XFRAC(1)»-1  >  XFRAC(2)=0.51 
XFRAC)3)s<*  J  XFRAC(<(  (=0.6475 
XFRACI5)=3  J  XFRACI6)=2.0 
XFRAC(7)=2  >  XFRAC(8)=I.4 


CROUP  4.  Y-OIRECnON  CRIO  SPECIFICATICN 
Ny*12 

yFRAC(ll«-3  t  yFRAC(2J«0.562 

yFRAC(S|s2  I  yFRACI4)«0.314 

yFRACI5)E3  t  yFRAC<4)>0.542 

yFRACI7)«4  I  yFRACIS)«0.7S 

GROUP  6.  Z'OIRECTZON  GRID  SPECIFICATION 


NZ<40 

ZFRAC(l)s-l 

1 

ZFRACi2)«1.0 

ZFRACI3I=1 

t 

ZFRAC(4)«I.7 

2FRACI5)>1 

> 

ZFRACI4)s2.0 

ZFRAC(7)3l 

1 

ZFRACI8)d.7 

ZFRACI9)s2 

1 

ZFRACd0)>3.0 

ZFRACdltsl 

> 

ZFRACdZ  t>1.5 

ZFRACI13)*! 

t 

ZFRACd4)>I.O 

ZFRACdS)*! 

» 

ZFRACd4)>I.7 

ZFRACd7)c2 

1 

ZFRACd6)«I.O 

ZFRACd9l«l 

1 

ZFRAC(  20)«1.7 

ZFRAC(21)s2 

) 

ZFRAC(22)d.O 

ZFRACC  23 (■3 

» 

ZFRACi 24 )>0. 564646 

ZFRAC( 25)^2 

1 

ZFRAC(24)«0.35 

ZFRAC(27)s7 

t 

2FRAC(28)s3.0 

ZFRAC(29)<1 

1 

ZFRAC(30)=2.0 

2FRAC(31)s4 

> 

ZFRAC(32)-0.7S 

ZFRAC(33)>4 

t 

ZFRACIS4)«0.575 

ZFRAC(3S)«4 

> 

ZFRAC(36)«0.5 

2FRAC(37)*1 

» 

ZFRAC(38)-3.0 

GROUP  4.  BOOY-FITTEO  COORDINATES  OR  GRID  DISTORTION 

GROUP  7.  VARIABLES  STORED*  SOLVED  t  NAMED 
S0LVE(P1,H1,U1,V1>HI) 

SOLUTNI  PI  * Y , Y , Y .N >N >N  ) 

ST0RE(T>1PI*RH01) 

GROUP  8.  TERMS  (IN  DIFFERENTIAL  EQUATIONS)  (  DEVICES 

GROUP  9.  PROPERTIES  OF  THE  MEDIUM 
TAMBS29S. 

TJETS1944. 

PJET=0. 

MJET*197.72 

PRESSO^IOISZS. 

RHOIA^O. 

GASC0N:286. 

RHOlBsl. /GASCON 

RHOJETst ( PRESSOtPJET )«RH01B )/TJET 

WJET=MJET/RHOJET 

TMPlA=nNY 

**»  TMPIB  •  1.  DIVIDED  BY  C>SUB-P)  USED  IN  TMPlsGRNOB 
CSUSP*1004. 

HAMB=TAMB»CSUBP 

TMPlEsl./CSUBP 

MJET=CSUBP«TJ£T 
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«•  NEXT  2  STATEHENTS  INVOKE  COHPRESSIBLE  FLQH  IN  CREXl 
TNP1«6RN02 
RH01*6RN05 
0RH10P*6RN05 
ENUU«l.E-05 
ENUr*200.«ENUL 
TURHOOCKEHOOLJ 
STORE! VIST) 

KE1NJ>0.5»I0.00S«HJET)*«2  |  EPINJ«0.09«KEINJ*«1.5/(0.0S«0.S) 
KEINA«0.5»(  .005*1.0 )M2  |  EPINA*0.09<HlEINAMl.S/i0.1) 

6RQUP  10.  INTER-PHASE 'TRANSFER  PROCESSES  AND  PROPERTIES 

GROUP  11.  INITIALIZATION  OF  VARIABLE  OR  POROSITY  FIELDS 
FIINITI HI )«HAMB 

FIIN1T!THP1)«TAMB  i  FIlNITIOENlJsl.O 

FIINIT(KE)«KEINA  >  FIINITI  EP  )>EP1NA 

FIINITI PI  )>0.0 

FIINITI HI  1*0.0 

FIINITI VI  )*0.0 

FIINITI  U1 1*0.0 

PATCH I HINIT , INI VAL ,1 >5  >1 >B  > 15  >30 , 1 , 1 ) 

INITI HINIT ,H1 ,0 . 0  >0 . 2«HJET I 
INITI HINIT>H1 >0.0  >CSUBP«0 . 2»TJET ) 

INITI NINIT  >KE  >0 . 0  >KE1NJ«0 . 2 ) 

INITI HINIT  >EP  >0 .0  >EP1NJ«0 . 2 ) 

•M  SEE  GROUP  13  FOR  POROSITY  FIELDS 

GROUP  12.  CONVECTION  AND  DIFFUSION  ADUUSTMENTS 

GROUP  13.  BOUNOARY  CONDITIONS  AND  SPECIAL  SOURCES 

«»0EFINE  SOURCES  AND  SII«S  FOR  JET  ENGINE** 
CONPORIO.O>CELL>1>1>4>5>12>14) 

PATCH! JETIN ,HIGH  >1>1>4>5>11>11>1>1) 

COVAL I JETIN , PI >  FIXFLU  t -MJET ) 

COVALI JEriN>Hl>ONLYHS>SAME I 
PATCH! JET0UT>L0W>1>1>4,5>1S>15>1>1) 

COVALI JET0UT>P1 , FIXFLU >MJET ) 

COVALI JETOUT,Hl>ONLYNS>HJET ) 

COVALI JET0UT>W1 ,ONLYMS>WJET ) 

COVALI JETOUT  >KE ,ONLyMS,KEINJ ) 

COVALI JETOUT,EP,ONLYMS>EPINJ ) 

*«EAST  BOUNOARY  TO  ATMOSPHERE** 

PATCHI  OUTSIDE, EAST>10>10,1>12>1>40,1>1) 

COVALI  OUTSIDE, PI  >0.1  >0.0) 

COVALI  OUTSIDE  >H1  ,FIXVAL>HA)1B  ) 

»*TOP  OF  HUSH  HOUSE  AND  TUBE  INCLUDING  TUBE  LIP** 

CONPORI 0.0, NORTH >1,9,12 >12, 1,19 1 

CONPORIO.O,NORTH,l,5,-8,-8,18,2r) 

COflPOR(0.0>NORTH,l,S,-2>-2>3S,3S) 

CO.NPOR  (  0 . 0  ,NORTH  ,  1 ,5  ,  -4  ,  -4 , 36  ,  34 ) 

CONPORI 0 . 0 ,NORTH ,1,5,-5,-5,37,371 
CONPORI 0.0>NORTH, 1,5, -6 ,-6,56, 38) 
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CONPORIO.O>NORTH.1.5>-7,-7*S9>S9l 


IMTOP  OF  MODEL  TO  ATMOSPHERE** 

PATCHC TOPOUr  ,NORTH,l >10 >12>12 > 20 »40 »1 ,1 ) 

COVAL I TOPOUr  >P 1 .0 . 1 >0 . 0 1 
COVALI  TOPOUr  ,H1,0NLYMS>HAIS  1 
PArCN<r0PP,N0RTH>10>10>12>12>1.19>l.IJ 
COVAL(TOPP>P1>0.1>0.0) 

COVAL I TOPP  >H1 >ONLVMS>HAMB ) 

mFLOOR  of  hush  HOUSEt  TUBE  AND  OROUNDtw 
CONPOR(0.0>SOUTH>1>10,1>-1>1,40) 

«*FORE  AND  AFT  BULKHEADS  IN  HUSH  HOUSE  AND  TUBE** 
CONPOR(0.0>LOH>1>9>1>12>-1>>1) 
CONPORCO.O>LOH,«>9,1>12>-20>-20) 
C0NPORI0.0,L0H.l>5>9>12.-20,-20) 
CONPOR(0.0>LOH>l,5>l>2>-55>-S5) 
C0NPOR(0.0>LOH>l>S,3>4>-36>-56) 

CONPORt 0.0, LOW, 1,5,5>5, -37,-37) 

CONPOR( 0.0, LOW, 1,S,A, 6, -30,-38) 
CONPOR(0.0,LOW,l,5,7,7,-39,-39) 

CONPOR( 0.0, HIGH, 1,5 ,8,8, -39,-39) 

*«VERY  FRONT  MODEL  STRIP  TO  ATMOSPHERE** 

PATCHI FRONT, LOW, 10, 10, 1,12, 1,1, 1,1) 

C0VAL(FR0NT,P1,0. 1,0.0) 

COVAL( FR0NT,H1,0NLYKS,HAM8) 

»»THE  VERY  ENO  OF  OUR  MODEL  TO  ATMOSPHERE** 

PATCH C EnOENO ,HI6H , 1 , 10 , 1 , 12 ,40 ,40 , 1 , 1 ) 

COVAL ( ENOENO ,P1,0. 1,0.0) 

COVALI ENOENO ,H1,0NLYMS,HAMB ) 

**OUTSIOE  EAST  )ULL  IN  HUSH  HOUSE** 
C0NP0RI0.0,EAST,-9,-9,9,12,l,I9) 
CONPORIO.O,EAST,-9,-9,l,8,l,l) 
CONPOR(0.0,EAST,-9,-9,l,8,3,3) 
CONPORIO.O,EAST,-9,-9,l,8,5,8) 
CONPORIO.O,EA3T,-9,-9,l,8,10,ll) 
CONPOR(0.0,EAST,-9,-9,l,8,13,14) 
CONPORIO.O,EAST,-9,-9,l,8,18,19) 

*«OUTSIOE  OF  TUBE  PLUS  TUBE  LIP  IN  HUSH  HOUSE** 
CONPORIO.O,EAST,-5,-5,l,8,18,39) 

CROUP  14.  OOM4STREAM  PRESSURE  FOR  PARAB«.TRUE. 

GROUP  15.  TERMINATION  OF  SHEEPS 
FSWEEPcl)  LSWEEP^IBO 

CROUP  16.  TER)1INATI0N  OF  ITERATIONS 

CROUP  17.  UNOER-RELAXATION  DEVICES 
REAL!  CELT, FACT) 
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DELT«SO . 0/( NZ«0 . ZkHJET  t 
FACT«10. 

FACT»2.S 

FACT«0.7S 

FACT>S.S 

FACT»S. 

FACT«O.SO 

FACT:0.2S 

RELAX(P1,LINRLX,0.3) 

RE  LAX( U1 , FALSOr , FACT«OELT ) 

RE  LAXI VI . FALSOT . F ACT«OE  LT ) 

RELAX! HI .FALSOT>FACT«OELT ) 

RELAX! HX,FALSOT>2.5«FACr»OELT) 

RELAX! KE , FALSOT, FACT»D£LT ) 

RELAX! EP. FALSOT, FACT*OELT) 

GROUP  18.  LIMITS  ON  VARIABLES  OR  INCREMENTS  TO  THEM 
VARMINIUDs'SOO. 

VARMAXIUI )slOOO. 

VARMZNIVl)*-300. 

VARMAXIVlIslOOO. 

VARMIN!Hl)s-100. 

VARMAX! HI  )>1500 . 

VARM1N!TMP1)=IS0 
VARMAX! TMPl »»TJET 
VARMIN!Hl)slOOO. 

VARMAX! HI  )sCSUBP*TJET 
VARMINIPl)3*3.5EtOS 
VARMAX! PI  S.SEtOS 
VARMlNIRHOlls.OS 
VARMAX! RH01)<5. 

VARMIN!KEJ=l.E-08 

VARMAX(KEI*1.5*K£JET 

VARMIN!EP)»l.E-08 

VARM1NIVIST)*0.1 

VARMAX!  EP  1*1.  S«(KEJET<Hil. 5) 

GROUP  19.  DATA  COMMUNICATED  BT  SATELLITE  TO  GROUND 

GROUP  20.  PRELIMINARY  PRIMT-OUT 

GROUP  21.  PRINT-OUT  OF  VARIABLES 

GROUP  22.  SPOT -VALUE  PRINT-OUT 

••MONITOR  THE  OUTPUT  AT  DESIGNATED  POINTS** 
IXMON^l  )  IYMON34  t  IZM0N=1A 
TSTSWP*!  I  LUPRS^A 
ABSIZ’l.O  >  QRSIZ^l.O 

GROUP  25.  FIELD  PRINT-OUT  AND  PLOT  CONTROL 
NUMCLS*10*KXPRIN 
ITABL*3  »  NPLT*1  >  NZPRIN=1 
IPLTLsLSWEEP 
NPRINTsS 
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GROUP  24.  DUMPS  FOR  RESTARTS 
SAVE ST 
RESTRT(AU) 

STOP 
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APPE34DXX  B 

SEL^EXn^ED  PHOEignrCS  OUTPUT 


CCCC  HHH 
CCCCCCCC  HHHHH 
CCCCCCC  HHHHHHHHHH 
CCCCCCC  HHHHHKHHHHHH 
CCCCCC  HHHHHHHHHHHHH 
CCCCCCC  HHHHHHHHHHHH 
CCCCCCC  HHHHHHHHHH 
CCCCCCCC  HHHHH 
CCCC  HHH 


PHOENICS  VERSXON  l.S,  03  SEPT  X986 
1C)  COPYRIGHT  1804 
C0)«:E))TRATZ0N  heat  atu  hohentum  lto 
ALL  rights  RESERVED. 

CHAM  LTD,  SAKERY  HOUSE,  40  HIGH  ST 
HIMOLEDON,  LONDON,  SN19  5AU 
TEL:  01-947-7451)  TELEX:  928517 
FACSIMILE:  01-879-3497 
THE  OPTION  LEVEL  IS  -18 


THIS  CODE  MAY  ONLY  BE  USED  UNDER  THE  TERMS  AND  CONDITIONS 
OF  A  LICENCE  AGREEMENT  HITH  CHAM  LTO. 

REPLICATION  OF  THIS  CODE  IS  PROHIBITED  UNLESS 


SPECIFICALLY  AUTHORISED  IN  WRITING  BY  CHAM  LTO. 


GREXl  OF  lS/07/86  HAS  BEEN  CALLED 

FORMATTED  SATLIT  DATA  READ  FROM  OFIO  FOR  IRUN*  1 

- - STORAGE  INFORMATION - - 


F  01MNSN=200000  0CCUPIE0»123707  ESTIMATED  MINIMUM  OIMNSN*  23247 
DEP  VRBL»  47200  OLD  VRSL*  0  30  COEFF»  25660  SO  DV0PS»  14400 


INITIAL  FIELDS  READ  FROM  0F09 
GROUP  1.  RUN  TITLE  I  NUMBER 

. . . . . -nnrnnnWMiiii»M»M»»M 


TEXTINAS  JACKSONVILLE  HUSH  CHAMBER 


IRUNN  X  X 

GROUP  2.  TRANSIENCE 
STEADY  X  T 

GROUP  3.  X-OIRE'CTION  GRID  SPACING 
CARTES  X  T 

NX  X  XO 

XULAST  X  1.000E*^00 

METHOD  OF  PAIRS  USED  FOR  GRID  SETTING 
XFRAC  (  X)  X  -I.000E,00  iXFRAC  (  2)  »  3.X00E-01 

XFRAC  (  3)  X  4.0002 too  )XFRAC  I  4)  x  5,475E-01 

XFRAC  I  S)  X  5.000E+00  )XFRAC  I  4)  «  2.0D0E*00 

XFRAC  1  7)  X  2.000E,00  )XFRaC  (  6)  ■  1.400E*00 

GROUP  4.  Y-OIRECTION  GRID  SPACING 
NY  X  X2 

YVLAST  X  X.OOOEtOO 

METHOD  Of  PAIRS  USED  FOR  GRID  SETTING 
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YFRAC 

( 

1) 

a 

-3.000E*00 

lYFRAC 

1 

Z) 

a 

S.A20E-01 

YFRAC 

( 

SI 

a 

2.aaoE*oo 

tYFRAC 

( 

4) 

a 

3.140E-01 

YFRAC 

( 

51 

a 

3.000E>00 

>YFRAC 

( 

A) 

a 

5.A20E-01 

YFRAC 

( 

7) 

a 

4.000E«00 

lYFRAC 

I 

8) 

a 

7.500e-01 

GROUP  5.  Z-OIRECTION  GRID  SPACING 
PARAB  «  F 
NZ  >40 


ZHLAST  >  l.OOOEtOO 

HETHOO  OF  PAIRS  USED  FOR  GRID  SETTING 


ZFRAC 

( 

11 

a 

-i.oaoE«oo 

IZFRAC 

( 

2) 

a 

1.000E»00 

ZFRAC 

( 

3) 

a 

1.000E*00 

iZFRAC 

I 

4) 

a 

1.700E»00 

ZFRAC 

( 

5) 

a 

1.000E>00 

iZFRAC 

( 

4) 

a 

2.000E«00 

ZFRAC 

( 

7) 

a 

1.000E«00 

)ZFRAC 

1 

B) 

a 

1.700E«00 

ZFRAC 

( 

91 

= 

2.000E«00 

tZFRAC 

( 

10) 

a 

i.OOOEtOO 

ZFRAC 

( 

11) 

s 

l.OOOE^OO 

)ZFRAC 

( 

12) 

a 

l.SOOE^OO 

ZFRAC 

( 

13) 

a 

1.000E«00 

IZFRAC 

( 

14) 

a 

l.OOOEtOO 

ZFRAC 

( 

IS) 

8 

l.OOOE+OO 

IZFRAC 

( 

16) 

a 

1.700E^00 

ZFRAC 

( 

17) 

8 

2.000E'»00 

tZFRAC 

( 

IB) 

a 

l.OOOE^CO 

ZFRAC 

( 

19) 

a 

l.OOOE^OO 

jZFRAC 

( 

20) 

a 

1.700E»00 

ZFRAC 

( 

21) 

a 

2.000E«00 

>2FRAC 

( 

22) 

a 

l.OOOEiOO 

ZFRAC 

( 

23) 

a 

S.OOOE^OO 

tZFRAC 

( 

24) 

a 

5.467E-01 

ZFRAC 

( 

25) 

a 

2.000E«00 

tZFRAC 

( 

24) 

a 

3.S00E-01 

ZFRAC 

( 

27) 

a 

7.000E«00 

tZFRAC 

( 

28) 

a 

S.OOOE^OO 

ZFRAC 

( 

29) 

a 

1.000£*00 

tZFRAC 

( 

30) 

a 

2 . OOOE  »00 

ZFRAC 

( 

31) 

a 

4.000E«00 

tZFRAC 

I 

32) 

a 

7.SOOE-01 

ZFRAC 

( 

33) 

a 

4.00QE*00 

tZFRAC 

( 

34) 

a 

S.7SOE-01 

ZFRAC 

( 

35) 

a 

4.000c*00 

tZFRAC 

I 

36) 

a 

5.000E-01 

ZFRAC 

( 

37) 

a 

l.OOOE+00 

tZFRAC 

( 

SB) 

a 

3. OOOE too 

CROUP  6.  BFC  AND  GEOMETRY  CHANGES 

GROUP  7.  VARIABLES:  STORED, SOLVED, NAMED 
ONEPHS  »  T 

NAME!  1)  »P1  JNAMEI  5)  *U1 
NAME!  5)  ■VI  )NAME(  7)  >N1 
NAME  (121  sKE  iNAMcdSI  «EP 
NAME (14 1  ■HI  )NAME(191  ■HPOR 
NAME(201  «NPCR  itUHSiZl)  *EPOR 
NAME (22)  «V?OR  iHAMEt231  *VIST 
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NAHE(24)  «RH01  )NAME(25)  >TMPX 
SOLUTNiPl  .Y,Y,Y,N>N>N) 

SOLUTN(Ul  ,Y.Y,N,N,N,N) 

SOLUTN(Vl  ,Y,Y,N,N,N>N) 

SOLinTKHl  ,Y.Y,N,N,N,N) 

SOLUT74(KE  ,Y«Y,N,N,N,N) 

SOLUrmEP  «Y,Y,N,N>N,N) 

SOLUTNIHI  ,Y>Y,N,N>N,NI 
SOLUTN(  HPOR,Y ,N,N,N,N>N ) 

SOLUTNi NPOR , Y ,N >N ,N ,N ,N ) 
S0LUTN(EP0R>Y,N,N>N,N,N1 
SOLUTNIVPOR,Y,N,N,N>N>N) 
SOLUTNIVIST,Y,N,N>N,N,N) 
S0LUTNIRH01,Y>N,N,N,N>N) 

SOLUTN( TMP1,Y ,N,N,N,N,N ) 

DENI  »  24 

VIST  «  23 

EPOR  3  21  (HPOR  «  19  »NPOR 

TEMPI  «  2S 


20  tVPOR 


GROUP  8.  TERMS  t  DEVICES 
TERMS  tPl  ,Y.Y,Y.N,Y,NI 
TERMS  tUl  ,Y,Y,Y,NiY.N) 

TERMS  (VI  ,Y,Y,Y,N,Y,N) 

TERMS  (HI  ,Y»Y,Y,N,Y,N) 

TERMS  (KE  ,N,Y,Y,N,Y,N1 
TERMS  (EP  ,N,Y,Y,N,Y,N) 

TERMS  (HI  ,Y»Y,Y,N,Y,N) 

OIFCUT  a  5.000E-01  jZDIFAC  «  X.OOOEtOO 
GALA  a  F  tAOOOIF  >  F 

NEWRHl  a  T 

NEWENT  a  T 


GROUP  9.  PROPERTIES 


RHOl 

a-1.0l6E*04 

jTMPX 

a-X.0X3E*04 

lELl 

a-1.01SE*04 

RHOIA 

a  O.OOOEtOO 

IRHOIB 

a  5.496E-03 

>RHOXC 

a  7.14$E-01 

PRESSO 

a  1.O13E+0S 

TMPIA 

a  l.pOOE-20 

iTMPie 

a  9.9o0E-04 

>TMP1C 

a  0.000E«00 

TEMPO 

a  O.OOOEtOO 

ELIA 

a  O.OCOE^OO 

,EL1B 

a  0.000E«00 

}EL1C 

a  O.OOOEtOO 
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ENUL 
EKUTA 
ORHIOP 
PRNOTLtUl 
PRNOTLINl 
PRNOTLIEP 
PRT  tUl 
PRT  IWl 
PRT  (EP 


s  l.OOOE-05  tENl/r  »1.014E«04 

■  0.000E«00  tENUTB  >  O.00OE*00  tENUTC  «  0.0006*00 


'1.016E«04 
) 

) 

) 


1.000E*00  >PRNDTUV1  )  •  1.000E«00 

1.000E400  tPRNDTUKE  )  *  1.000E»00 

l.OOOE^OO  iPRNDTLCHl  )  >  1.000E«00 

)  s  1.000E«00  »PRT  fVl  )  B  l.OOOE^OO 

)  *  1.000E*00  )PRT  IKE  )  ■  X.OOOE^OO 

)  *  1.314Et00  tPRT  (HI  J  «  1.000E«00 


GROUP  10. INTER-PHASE  TRANSFER  PROCESSES 


CROUP  11. INITIALIZE  VAR/POROSITY  FIELDS 
FIINITCPI  )  »  -1.02ZE*04  »FIINIT(U1  )  «  -1.022E«04 

FIINITIVI  )  *  -1.022E«04  IFIINITIHI  )  «  -1.022E«04 

FIINITIKE  }  «  -1.022£«04  >FIINIT1EP  )  *  -1.022£*04 

FIINITIHl  )  «  -1.022E+04  jFIINITI HPOR )  a  -1.022E+04 
FIINIT(NPOR)  a  -l,022£+04  IFIINITI EPOR )  a  -1.022E+04 
FIINIT(VPOR)  a  -1,022E*04  iFIINITCVIST)  a  -1.022E+04 
FIINITIRHOl)  a  -X.022E+04  >FIINITITHP1)  «  -1.022E+04 


PATCHI SOUTH! 

INITISOUTHl 

INITISOUTHl 

PATCH! LOHl 
IN1T(L0H1 
INITI LOWl 

PATCH! CELLl 

INIT!CELL1 

INIT!CELL1 

INITICELLI 

INIT!CEIL1 

INITICELLI 

INITICELLI 

INITICELLI 

INITICELLI 

INITICELLI 

INITICELLI 


>INIVAL>  1,  1>  3»  3, 


rVl  , 

0.000E>00, 

0.000E«00) 

iNPOR, 

0.00OE«0O, 

O.OOOE+OO) 

,INIVAL,  1, 

1,  4,  5 

,W1  , 

0.000E«^00, 

O.OOOEtOO) 

.HPOR, 

O.OOOEtOO, 

O.OOOEtOO) 

,INIVAL,  1, 

1,  4,  5 

PI  , 

O.OOOE^OO, 

O.OOOETOO  ) 

.U1  , 

0.000E«00, 

O.OOOEtOO) 

VI  , 

0.00OE*OO, 

O.OOOEtOO  1 

.W1  , 

0.OO0E*OO, 

O.OOOEtOO  ) 

iKE  , 

O.OOOEtOO, 

O.OOOEtOO  ) 

EP  , 

O.OOOE^OO, 

O.OOOEtOO) 

.HI  , 

O.OOOE^OO, 

O.OOOEtOO  ) 

iHPOR, 

O.OOOE*00, 

O.OOOEtOO) 

iNPOR, 

O.OOOE+OO, 

O.OOOEtOO ) 

EPOR, 

O.OOOE»00, 

O.OOOEtOO) 

K>  14» 

11>  11, 

IZ,  14, 


1, 


1, 


1, 


1) 


1) 


1) 
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5 


5 


1 


1) 


PATCH(EASrzi  ,INZVAL,  9,  9,  1,  B, 

INlT(EASr21  ,U1  f  O.OOOE^OO.  O.OOOEtOO) 

ZNITIEAST21  ,EPOR,  O.OOOE«00>  O.OOOE^OO) 

PATCH!  EAST22  .INIVAL,  9,  9,  1,  8>  S,  8>  1,  1) 

INIT(EAST22  >U1  ,  O.OOOE«00>  O.OOOE^OOJ 

INIT(EAST22  >EPOR>  0.000E«00>  O.OOOE*00) 

PATCH!  EAST23  ,INIVAL.  9,  9,  1,  6,  10,  11,  1,  1) 

1NIT!EAST23  ,U1  ,  O.OOOE^OO,  O.OOOE^OO) 

INIT!EAST23  >EPOR>  O.OOOEtOO,  O.OOOE^OO) 

PATCH!EAST24  ,IN1VAL>  9,  9,  1,  6,  13,  14,  1,  1) 

INIT!EAST24  lUl  ,  O.OOOE^OO,  Q.OOOttOO) 

INIT(EASr24  ,EPOR,  0.000£*00,  0.000E*00) 

PATCH! EAST25  ,IN1VAL,  9,  9,  1,  8,  18,  19,  1,  1) 

INITIEAST25  ,U1  ,  0.000E«00,  O.OOOE«00) 

INIT!EAST25  ,EPOR,  O.OOOEtOO,  O.OOOE'fOOJ 

PATCH! EAST26  ,INIVAL,  S,  S,  1,  8,  18,  39, 

INIT!EAST26  ,U1  ,  O.OOOE*00,  O.OOOE^OO) 

1N1T!EAST26  ,EPOR,  O.OOOEtOO,  O.OOOEtOO) 

INIAOO  «  F 
LUFI  «  9 

NAMFI  *OF09 

GROUP  13.  BOUNDARY  i  SPECIAL  SOURCES 

PATCH!  KESOURCE,PHASEri,  1,  10,  1,  12,  1,  40,  1,  1) 

COVALIKESOURCEjKE  , -I . 015E»04,-1 .OlSE+04 ) 

C0VAL1KES0URCE,EP  ,-1.01SE*04,-1.01SE+04 ) 

PATCH! JETIN  .HIGH  ,  1,  1,  4,  5,  11,  11,  1,  1) 

COVALIJETIN  ,Pi  ,  1.000E-10,-1.977Et02) 

COVALI  JETIN  ,H1  ,  O.OO0E«0O,-1.026EtO4 ) 

PATCHIJETOUT  ,LOW  ,  1,  1,  4,  5,  15,  15,  X,  1) 

COVAL!JETOUT  ,P1  ,  I.OOOE-IO,  1.977E*0Z) 
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C0VALIHWALL12 

>U1 

,-1.013E«^04.  0.000E«00) 

C0VALIHMALL12 

,V1 

.-1.013E+04,  O.OOOEfOO) 

C0VALIHHALL12 

*KE 

,-1.013£«04,-1.013E+04) 

COVALIHHALLia 

»EP 

.•1.013E«04.-1.013E«04) 

PATCH(LMALL12 

,LHALL  ,  1,  S,  12, 

20. 

20, 

1. 

1) 

COVAD  LHALL12 

>U1 

.-1.013E«04.  0.000E«00) 

COVAD  I.HALL12 

,V1 

.-1.013E»04.  0.000E«00) 

COVAD  LHALL12 

,KE 

.-1.013E404.-1.013E'»04) 

COVAL( LHALL12 

,EP 

.-1.013E«04.-1.013Et04) 

PATCH(HHALL13 

,HHALL  ,  1.  S,  1,  2, 

34. 

34, 

1. 

1) 

COVAD  HHALL13 

,U1 

.-1.013E*04.  0.000E«00) 

C0VAL(HHALL13 

,V1 

.-1.013E»04.  0.000E«00) 

COVAD  HMALL13 

,KE 

.-1.013Et04.-1.013E*04) 

COVAD  HHALL13 

,EP 

.-1.013E*04.-1.0i3E04) 

PATCH( LWALL13 

,LWALL  ,  1,  5,  1,  2, 

35. 

35, 

1. 

1) 

COVAD  LHALL13 

,U1 

.-1.013E»04,  0.000E400) 

COVAD  LHALL13 

.VI 

.-1.013E404.  O.OOOE+OO) 

COVAD  LMALL13 

,KE 

.-1.013E404,-1.013E404) 

COVAD  LWALL13 

,EP 

.-1.013E«04.-1.013Et04) 

PATCH(HHALL14 

»HWALL  p  X*  5» 

35. 

35. 

1, 

1) 

C0VAL(HHAD.14 

,U1 

.-1.013E404,  O.OOOE400) 

COVAD  HI-iALLU 

.VI 

.-1.013E404.  0.000E400) 

COVAD  H»^(ALL14 

,KE 

.-1.015E404.-1.013E 404 ) 

COVAD  HWALL14 

.EP 

.-1.013E404.-1.013E404 ) 

PATCH(LWALL14 

.IHALL  .  1.  5.  3.  4. 

36. 

36, 

1. 

11 

COVAD  LHALL14 

.01 

.-1.013E404.  0.000E400) 

COVAD  LHALL14 

.VI 

.-1.013E404.  O.OO0E4OO) 

COVAD  LNALL14 

.KE 

,-1.013E404.-1.015E404) 

COVAD  LHALL14 

.EP 

.-1.013E404.-1.013E404) 

PATCH! HHALL15 

fHWALL  f  Xp  5p  5p  5p 

36, 

36, 

1. 

1) 

COVAD  HHALL15 

.U1 

.-1.013E404,  0.000E400) 

COVAD  HHALL15 

.VI 

.-1.013E404,  O.OOOE400) 

COVAD  HKALL15 

,KE 

.-1.013E404.-1.01IE404) 

COVAD  KHALLIS 

.EP 

.-1.013E404.-1.C1$E404) 

I 


I 
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PATCH! LNALL15 

,LNALL  «  1.  5,  5.  S, 

57. 

57, 

1, 

1) 

COVALI LHALL15 

»U1 

>-1.015E*04,  O.OOOEOO) 

COVAL! LHALLI5 

tVl 

.-1.01SE«08>  O.OOOEOO) 

COVALI LHALL15 

>KE 

>-1.013E«04.-1.0XSE«04 ) 

COVALI  LHALL15 

.£P 

, -1 . 012E *04 , -1 . 015E *04 ) 

PATCHIHWALL16 

>HKALL  ,  1,  5,  8,  8, 

57, 

57, 

1, 

1) 

COVALI HMALLIA 

tUl 

.-1.0UE404,  0.000E*00) 

COVALI HHALL16 

.VI 

.-1.01SE«04,  O.OOOE^OO) 

COVALI HMALL16 

.KE 

.-1.013E«04,-1.013E«04) 

COVALI HHALL16 

.EP 

,-1.015E*04,-1.015E*04J 

PATCHI LHALL16 

.LHALL  ,  1,  S,  8.  8. 

58, 

58, 

X, 

X) 

COVALI  LHALL16 

.U1 

.-1.0I3E404,  O.OOOE400) 

COVALI LHALL16 

.VI 

.•1.013E404.  0.000E«00) 

COVALI LWALL16 

.KE 

.  -1 . 013E  *04 . -1 , 013E  *04 ) 

COVALI LHALL16 

.EP 

.-1.013E404.-X.013E 404) 

PATCHI HWALL17 

.HHALL  >  1.  5.  7.  7. 

58. 

58, 

1, 

X) 

COVALIHWALL17 

.U1 

.•X.0X3E404.  O.OOOE400) 

COVALI HHALL17 

.VI 

.-1.013E404.  O.OOOE40O) 

COVALI HHALLl 7 

.KE 

.-X.0X3E404.-1.0X3E404) 

COVALI HHALL17 

.EP 

, -1 . 0X3E  404 , -1 . 013E  404 ) 

PATCHI LMALL17 

.LHALL  .  1.  5.  7.  7. 

39. 

59, 

1, 

1) 

COVALI LHALL17 

.U1 

.•X.0X3E404,  0.000E400I 

COVALI LHALL17 

.VI 

,-X.013E404.  O.OOOF40O) 

COVALI LHALL17 

.KE 

.-X.013E404.-1.013E404) 

COVALI LHALL17 

.EP 

>-X.0X3E404,-X.0X3E404) 

PATCHI HHALL18 

pHHALL  9  Xp  5p  Sp  8p 

59. 

59, 

X, 

X) 

COVALI HHALL18 

,U1 

.-X.013E404.  0.000E400) 

COVALI HHALL18 

.VI 

.-X.013E404,  O.000E400) 

COVALI KHALL18 

.KE 

,-1.0X3E404.-1.013E404) 

COVALI HWALLle 

.EP 

,-1.0I3E404.-X.0X3E404) 

PATCHI LHALLia 

pLWALt  p  Xp  5p  8p  Sf 

40. 

40, 

X, 

XI 

COVALI LHALL18 

,U1 

.-1.0X3E404.  O.O00E400I 

COVALI LHALL18 

.VI 

.-X.0X3E404.  O.O00E400) 
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COVALI LHALL18 

»KE 

.•1.013E*04.-1.01SE+04) 

COVAL( LHALL18 

.EP 

.-1.013E«04.-1.013E»04) 

PATCHr  FRONT 

>LOH 

.  10.  10.  1,  IZ. 

1. 

1. 

1. 

1) 

COVALI FRONT 

,P1 

>  l.OOOE-01.  O.OOOE^OO) 

COVAL(  FRONT 

>H1 

,  0.000£«00.  Z.96ZE«05) 

PATCH!  ENDENO 

.HIGH 

.  1.  10.  1,  IZ. 

40. 

40. 

1. 

11 

COVAKENDENO 

.PI 

.  l.OOOE-01.  O.OOOE.OO) 

COVALC  ENOEtJO 

.HI 

.  0.000E«00.  Z.96ZE4^05) 

PATCH! EMALL19 

.EHALL  .  9.  9.  9.  IZ, 

1. 

19. 

1, 

1) 

COVAL! EWALL19 

.VI 

.-1.013£«04,  O.OOOE+00) 

COVAL!  EHALL19 

.HI 

.-1.013E«04.  O.OOOE.OO) 

COVAL!  EHALL19 

.KE 

.-1.013E.04.-1.013E.04 ) 

COVAL! EHALL19 

.EP 

.•1.013E»04.-1.013E»04) 

PATCH!  FMALL19 

.('MALL  .  10.  10.  9.  IZ, 

1. 

19. 

1. 

1) 

COVAL! HMALL19 

.VI 

.>1.013E*04,  O.OOOEOO) 

COVAL!  WWALL19 

>H1 

.•1.013E*04,  O.OOOE.OC) 

COVALI  HHALL19 

>KE 

,-1.013E«04.>1.013E.04l 

COVAL!  WULL19 

.EP 

. -1 . 013E  «  04 . -1 . 013E  t04 ) 

PATCH! EHALLZO 

.EHALL  .  9.  9.  1.  8. 

1, 

1. 

1. 

1) 

COVAL!  EHALLZO 

.VI 

.-1.013E*04.  O.OOOEfOO) 

COVAL! EHALLZO 

.HI 

.•1.013E«^04.  O.OOOEtOO) 

COVAL! EHALLZO 

.KE 

.-1.013E*04.>1.013Ef04) 

COVAL! EHALLZO 

.EP 

.-1.013E«04,-1.013Et04) 

PATCH!  HHALLZO 

.HHALL  .  10,  10,  1,  8, 

1. 

1, 

1, 

1) 

COVAL!  HHALLZO 

.VI 

.-1.013Et04,  O.OOOEtOO) 

COVALIHHALLZO 

.HI 

.-1.013Et04.  O.OOOEtOO) 

COVAL! HHALLZO 

.KE 

.-1.013Et04,-1.013Et04 ) 

COVALIHHALLZO 

.EP 

. -1 . 013E  t04 . -1 . 013E  t04 ) 

PATCH! EHALLZl 

lEWALL  f  9p  It  8t 

3, 

3. 

1. 

1) 

COVAL! EHALLZl 

,V1 

,-1.013Et04,  O.OOOEtOO) 

COVAL!  EHALLZl 

.HI 

.-1.013Et04,  O.OOOEtOO) 

COVAL! EHALLZl 

,KE 

.-1.013Et04,-1.013Et04 1 

COVAL!  EHALLZl 

.EP 

.-1.013Et04,-1.013Et04) 
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PATCH) HMALL21 

,HWALL  >  10,  10,  1,  a. 

3, 

3, 

1, 

1) 

COVAD  l«IALL21 

,V1 

,-1.013E«04,  O.OOOE^'OO) 

COVAD  l>HALL21 

,M1 

,-1.013E«04,  0.O0OE«^OO) 

COVAD  HHAD.21 

,KE 

,-1.013E«04,-1.013E«04 ) 

COVAD  NWALL21 

,EP 

,-1.013E404,-1.013E«04) 

PATCH) EHALL22 

,E»)AU  ,  9,  9,  1,  a. 

5, 

8, 

1, 

1) 

COVAD  EliALL22 

,V1 

,-1.013E404,  O.OOOE^OO) 

COVAD  EHALL22 

,W1 

,-1.013E»04,  O.OOOE«00) 

COVAD  EHALL22 

,KE 

,-1.013E*04,>1.013E04 ) 

COVAD  EHALL22 

,EP 

,*1.013E*04,-1.013Et04) 

PATCH) HHALL22 

,HHALL  ,  10,  10,  1,  8, 

5, 

8, 

1, 

1) 

COVAD  l«ULL22 

,V1 

,-1.013E*04,  O.000E«00) 

COVAD  HHALL22 

,H1 

,>1.013E«04,  O.OOOE«^O0> 

COVAD  HMALL22 

,KE 

,>l.ai3E«04,-1.013E«04) 

COVADKHALL22 

,EP 

,-1.013E*04,*1.013Et04) 

PATCH) EWALU23 

, EMAIL  ,  9,  9,  1,  8, 

10, 

11, 

1, 

1) 

COVAD  EHALL23 

,V1 

,>1.013Et04,  0.000E«00) 

COVAD  ENALL23 

,W1 

,-1.013E*04>  O.OOOE^OOJ 

COVAD  EHALL23 

,KE 

•-1.013E«04,-1.013Et04) 

COVAD  EHAD.23 

,EP 

, *1 . 013E«04, -1 . 013E  t04 ) 

PATCH) HHALL23 

,MHALL  ,  10,  10,  1,  8, 

10, 

11, 

1, 

1) 

COVAD  HHALL23 

,V1 

,*1.013E4’04,  O.OOOEiOO) 

COVAl.)WWALL23 

,HX 

,-1.013E«04,  O.OOOEtOO) 

COVAD  HHALL23 

,KE 

,-1.013Et04,-1.013E*04) 

COVAD  K-)ALL23 

,EP 

,>1.013E*04,-1.013EtC4 ) 

PATCH) EHALL24 

, EMAIL  ,  9,  9,  1,  8, 

13, 

14, 

1, 

1) 

COVAD  EHALL24 

,V1 

,*1.013E»04,  O.OOOEOOl 

COVAD  EWALL  24 

,M1 

,>1.013E*04,  O.OOOEfOO) 

COVAD  EHALL24 

,KE 

,'1.0i3E+04,-1.013E»04) 

COVAD  EHALL24 

,EP 

,-1.013E»04,-1.013Et04) 

PATCH) KHALL24 

,WMALL  ,  10,  10,  1,  8, 

13, 

14, 

1, 

1) 

COVAD  WMALL24 

,V1 

,-1.013Et04,  O.OOOEtOO) 

COVAD  )<XALL24 

,W1 

,-l.ai3Et04,  O.OOOEtOO) 
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C0VAL(MHALL24  ,KE  ,-1.015E*04,-X.01SE»04) 

C0VAL(WWALL24  ,EP  ,-1.015E*04,-1.015E*04 ) 

PATCHIEHALL2S  ,EHALL  ,  4,  4,  1,  8>  18,  14,  1,  1) 

COVAUt EHALL2S  ,V1  ,-1.013E*04,  O.OOOE^OO) 

COVAL(EWALL2S  ,N1  ,-1.013E«04,  O.OOOE+00) 

COVAL( EWALL25  ,KE  ,-1.013E«04,-1.013E«04 ) 

COVAL(EHALL2S  ,EP  ,-1.013Ei04,-1.013E*04 ) 

PATCH(WWALL25  ,HWALL  ,  10,  10,  1,  8.  18,  19,  1,  1) 

C0VALJ»*<ALt2S  ,V1  ,-1.013E«04.  O.0O0B*00) 

C0VAL(HWALL25  ,H1  ,-1.013E*04,  O.000E«00) 

C0VAL(HHALL2S  ,KE  ,-1.013E*04,-1.013E*04) 

C0VAL(HHALL2S  ,EP  ,-1.013E»04,-1.015E*04 ) 

PATCHt EWALL26  .EWALL  ,  5,  5,  1,  8,  18,  54,  1,  1) 

COVAL(EWALL26  ,V1  ,>1.013E'»04,  O.OOOE«00) 

COVAL(EHALL26  ,H1  •-1.013E*04,  O.OO0E«OO) 

C0VAL(EWAUL26  ,KE  ,-1.013E*04,-1.015Ei04 ) 

COVAL(EHALL24  ,EP  ,-1.013£»04,-1.013E*04) 


PArCH(WWAU24  ,HKAU  ,  4,  4,  1,  8,  18,  34,  1,  1) 

C0VAL(HHALL26  ,V1  ,-1.013E»04,  O.OOOE+00) 

COVAL(MWALL24  ,W1  ,-1.013E*04,  O.OOOEtOO) 

C0VAL(HWALt.24  ,KE  ,-1.013E*04,-1.01$E»04) 

COVALIHWALL24  ,EP  ,-1.013E^04,-1.013Ef04 ) 

XCYCLE  ■  F 
KELIN  <  0 


GROUP  14.  DOWNSTREAM  PRESSURE  FOR  PARAB 


GROUP  IS.  TERMINATE  SWEEPS 


LSWEEP  s  180  )ISWC1  «  1 

LITHYD  »  1  jLITFLX  «  1  »LITC  »  1  jITHCl  « 

RESREFIPI  I  *  l.OOOE-06  iRESREFlUl  i  -  l.OOOE-04 

RESREF(V1  )  «  1.000E-06  (RESREFIHl  )  «  l.OOOE-04 

RESREFIKE  )  *  l.OOOE-06  (RESREFIEP  )  *  l.OOOE-06 

RESREFOIl  )  s  l.OOOE-06 


GROUP  16.  T£R;iINATE  ITERATIONS 
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LITER 

(PI 

) 

S 

20  }LITER  (U1  )  « 

20 

LITER 

(VI 

) 

s 

20  >LITER  (Ml  )  > 

20 

LITER 

(KE 

) 

s 

20  tLITER  (EP  )  « 

20 

LITER 

(HI 

) 

s 

20 

ENOIT 

(PI 

) 

8 

l.OOOE-06  tENOIT  (U1 

)  » 

l.OOOE-06 

ENOIT 

(VI 

) 

8 

l.OOOE-06  tENOIT  (Ml 

J  » 

l.OOOE-06 

ENOIT 

(KE 

) 

8 

l.OOOE-06  tENOIT  (EP 

)  » 

l.OOOE-06 

ENOIT 

W  W  W  W  fc. 

(HI 

) 

ftIH 

8 

»«« 

1.000E-06 

CROUP  17.  RELAXATION 


RELAX(P1  ,LINRLX,  3.000E-01) 
RELAXfUl  ,FALS0T,  4.321E-0S) 
RELAX(V1  ,FALS0r,  4.321E-03) 
RELAX(H1  .FALSOT,  4.321E-03) 
RELAX(KE  , FALSOT,  4.321E-03) 
RELAX(EP  , FALSOT,  4.321E-03) 
RELAX(H1  , FALSOT,  1.080E>02) 
RELAX( VIST, FAt SOT,  1.000E*00» 
RELAX(RK01, FALSOT,  l.OOOE^OOl 
RE LAX(TMP1, FALSOT,  1.000E*00) 


GROUP  18.  LIMITS 

VARMAXIPl  i  •  3.500E«0S  iVARMINIPl  )  *-3.500E«05 
VARMAXIUl  )  *  1.000E*0S  »VARMIN(U1  )  «-3.000E*02 
VARMAXIVl  )  =  l.OOOE+03  »VARMIN(V1  )  s-3.000E*02 
VARMAXCHl  )  :  1.500E«03  iVARMINIHl  }  >-1.000E«02 
VARMAXIKE  )  =  l.OOOEtlO  jVARMINIKE  J  »  l.OOOE-10 
VARMAX(EP  )  =  l.OOOEilO  iVARMINtEP  )  «  l.OOOE-10 
VARMAX(H1  )  =  1.952Et04  tVARMINIHl  )  *  1.000E«03 
VARMAXfHPORI  <  l.OOOEtlO  }VARMIN( HPOR )  >-1.000E*10 
VARMAX(NPORI  s  I.OOOE^IO  tVARMIN( NPOR )  «-1.000E«10 
VARMAXIEPORJ  =  l.OOOEtlO  iVARMINI EPOR )  *-1.000£*10 
VARMAX(VPOR)  »  l.OOOEtlO  iVARMINt VPOR )  «-1.000E+10 
VARMAX(VIST)  s  l.OOOE+10  jVARMINI VIST )  «  l.OOOE-20 
VARMAXIRHOl)  s  S.OOOE^OO  iVARMIHI RHOl)  =  5.000E-02 
VARMAXITMPIJ  >  1.944E«03  )VARMIN( TMPl )  s  2.500Et02 


GROUP  19.  EARTH  CALLS  TO  GROLPiO  STATION 
USEGRO  >  T  tUSEGRX  s  T 
NAMGRO  sNONE 
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CROUP  20.  PRELIMINARY  PRINTOUT 
ECHO  >  T 

CROUP  21.  PRINT -OUT  OF  VARIABLES 
INI FLO  ■  F  jSUBNCR  «  F 

OUTPUTtPl  ,Y,N,N,Y,Y,Y» 

OUTPUTIUI  ,Y,N,N,Y,Y,Y) 

0UTPUT(V1  ,Y,N,N,Y,Y,Y) 

0UTPUT(W1  ,Y,N,N,Y,Y,YI 
OUTPUTCKE  ,Y,N,N,Y,Y,YI 
OUTPUTIEP  ,Y,N,N,Y,Y,Y) 

OUTPUT*  HI  ,Y,N,N,Y,Y,Y) 

OUTPUT* HPOR>N,N,N.N,N,NI 
OUTP*JT  ( NPOR  ,N  ,N  ,N  ,N,N  >N ) 

OUTPUT* EPOR>N,N,N,N,N,N) 

OUTPUT* VPOR>N,N,N,N,N>N) 

OUTPUT* VIST, Y,N,N,N,N,N) 

OUTPUT*  RH01,Y,N,N,N,N,N) 
0UTP*JT*TMP1,Y,N,N,N,N,NI 

GROUP  22.  MONITOR  PRINT -OUT 
IXMON  s  1  jIYMON  »  ^  fIZMON  * 

NPRMON  *  10000  JNPRMNT  *  10000  jTSTSWP  » 
HIGHLO  «  F 

GROUP  23. FIELD  PRINT -OUT  t  PLOT  CONTROL 


NPRINT 

m 

10000 

iNUMCLS 

8 

10 

NXPRIN 

s 

1 

lIXPRF 

8 

1 

JIXPRL 

8 

10 

NYPRIN 

m 

1 

JIYPRF 

8 

1 

JIYPRL 

8 

12 

NZPRIN 

X 

1 

lIZPRF 

8 

1 

lIZPRL 

8 

40 

XZPR 

8 

F 

lYZPR 

8 

F 

IPLTF 

8 

1 

JIPLTL 

8 

ISO 

»NPLT 

8 

1 

ITABL 

8 

3 

tIPROF 

8 

1 

ABSIZ 

s 

l.OOOEtOO  lORSIZ 

»  1. 

OOOEtOO 

NT2PRF 

8 

1 

tNCOLPF 

8 

50 

ICHR 

8 

2 

iNCOLCO 

8 

<»5 

(NROHCO 

8 

20 

NO  PATCHES  YET  USED  FOR  THIS  GROUP 
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GROUP  24.  DUMPS  FOR  RESTARTS 
SAVE  «  T  lAUTOPS  «  F  INOWIPE  ■  F 
NSAVE  ■0F09 
LUSAVE  «  9 

4HHf  GRID-GEOMETRV  INFORMATION  ••• 


X-COOROINATES 

OF  THE  CELL 

CENTRES 

l.SSOE-01 

S.837E-01 

1.131E«00 

1.679E«^00 

S.SOOE^OO 

S.SOOEtOO 

7.500EO0 

9.200E«00 

Y-COOROINATES 

OF  THE  CELL 

CENTRES 

2.81.0E-01 

8.430E-01 

1.405E*00 

1.843E«00 

2.59SE*00 

3.1S7E«00 

5.719E+00 

4.375EOO 

5.875E*00 

6.62SE*00 

Z-COOROINATES 

OF  THE  CELL 

CENTRES 

5.000E-01 

I.8S0E»00 

S.700E*00 

5.550E*^00 

1.090E401 

1.315E401 

1.440E*0I 

1.575E*01 

1.810£>01 

1.94Sc*01 

2.060E*04. 

2.180E*01 

2.31SE«01 

2.372E*01 

2.417E*01 

2.4S2r  "L 

2.920E*01 

3.220Et01 

3.S20E«^01 

3.820E«0I 

4.420E+01 

4.670E«01 

4.807E»01 

4.8d2Ei01 

S.022Et01 

5.099E+01 

S.156E*01 

5.214E«01 

S.325E>01 

5.37£Et01 

5.42SE»01 

5.475E»01 

-  INTEGRATION  OF  EQUATIONS  BEGINS  — 

_ 

TIME  STEP  *  1 

SHEEP  « 

2 

TOTAL 

RESIDUAL/I 

I.OOOE-061 

FOR 

PI 

IS 

8.007E*05 

TOTAL 

RESIDUAL/I 

l.OOOE-06) 

FOR 

U1 

IS 

9.092Et06 

TOTAL 

RESI0UAL/( 

1.000E-a6) 

FOR 

VI 

IS 

1.282E+07 

TOTAL 

RESIDUAL/I 

l.OOOE-06) 

FOR 

HI 

IS 

2.37fcE*07 

TOTAL 

RESIOUAL/I 

l.OOOE-06) 

FOR 

KE 

IS 

1.210E«09 

TOTAL 

RESIDUAL/I 

l.OOOE-06) 

FOR 

EP 

IS 

1.49bE+ll 

TOTAL 

RESIDUAL/t 

l.OOOE-06) 

FOR 

HI 

IS 

5.642E«10 

TIME  STEP  =  1 

SWEEP  » 

3 

TOTAL 

RESIDUAL/( 

l.OOOE-06) 

FOR 

PI 

IS 

3.7S3Ef05 

TOTAL 

RESIOUAL/I 

l.OOOE-06) 

FOR 

U1 

IS 

1.0£0E«07 

TOTAL 

RESIDUAL/I 

l.OOOE-06) 

FCR 

VI 

IS 

1.474E*07 

TOTAL 

RESIOUAL/I 

l.OOOE-06) 

FOR 

HI 

IS 

2.057E+07 

TOTAL 

RESIOUAL/I 

l.OCOE-04) 

FOR 

KE 

IS 

1.424E*09 

2.226Et00 

1.060E«01 

2.157E«00 

5.125E*00 


7.900E»00 
1.710Et01 
2.256E4^01 
2.620E4^01 
A.lZOEtOl 
4.9S7E  *01 
5.271E«01 
5.650E4^0I 
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TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
HUE  STEP  ■  1 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TIME  STEP  »  1 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TIME  STEP  «  1 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TIME  STEP  s  1 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/I 
TOTAL  RESIDUAL/I 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/I 
TOTAL  RESIDUAL/! 
TIME  STEP  »  1 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 
TOTAL  RESIDUAL/! 


1.000E-06) 

FOR 

EP 

l.OOOE-Ob) 

FOR 

HI 

SWEEP  « 

174 

l.OOOE-06) 

FOR 

PI 

l.OOOE-06) 

FOR 

U1 

l.OOOE-06) 

FOR 

VI 

l.OOOE-061 

FOR 

HI 

l.OOOE-Ob) 

FOR 

KE 

l.OOOE'06) 

FOR 

EP 

I.OOOE-06) 

FOR 

HI 

SWEEP  ■ 

17S 

l.OOOE-06) 

FOR 

PI 

I.OOOE-Ob) 

FOR 

U1 

l.OOOE-06) 

FOR 

VI 

l.OOOE-061 

FOR 

HI 

l.OOOE-06) 

FOR 

KE 

l.OOOE-Ob) 

FOR 

EP 

I.OOOE-06) 

FOR 

HI 

SHEEP  « 

17b 

l.OOOE-Ob) 

FOR 

PI 

1,OOOE-D6) 

FOR 

U1 

l.OOOE-Ob) 

FOR 

VI 

l.OOOE-06) 

FOR 

HI 

I.DOOE-06) 

FOR 

KE 

I.OOOE-Ob) 

FOR 

EP 

I.000£-06i 

FOR 

HI 

SWEEP  » 

177 

1.  CODE -Ob) 

FOR 

PI 

l.OOOE-Ob) 

FOR 

U1 

l.OOOE-Ob) 

FOR 

VI 

l.OOOE-Ob) 

FOR 

HI 

l.OOOE-Ob) 

FOR 

KE 

l.OOOE-Ob) 

FOR 

EP 

l.OOOE-Ob) 

FOR 

HI 

SHEEP  ■ 

178 

l.OOOE-Ob) 

FOR 

PI 

l.OOOE-Ob) 

FOR 

U1 

l.OOOE-Ob) 

FOR 

VI 

l.OOOE-Ob) 

FOR 

HI 

l.OOOE-Ob) 

FOR 

KE 

ZS  2.S41E«10 
IS  2.A20E*1D 

ZS  5.251E«04 
ZS  A.D2IE«D6 
ZS  6.189E»06 
ZS  6.271E«06 
IS  2.926Et0a 
IS  2.45£E»10 
IS  2.A19E«^10 

IS  S.208Et04 
IS  5.996EiD6 
IS  A.iaiEtDb 
ZS  8.329E406 
IS  2.91£E«08 
IS  2.393EtlO 
ZS  2.62aE»ID 

ZS  3.193Ei04 
ZS  6.D71EiD6 
IS  6.125E«06 
IS  8.259E«06 
IS  2.909E«08 
IS  2.4S5EilD 
IS  Z.6Z8£*lO 

IS  3.235Ei04 
ZS  6.038E»06 
ZS  6.I28EtD6 
IS  8.291E»06 
IS  2.9S4E«08 
IS  2.4b6EtlO 
IS  Z.622EtI0 

IS  S.272E*04 
IS  b.095Et06 
IS  6.115Et06 
IS  D.309Et06 
IS  S.DlDEfOa 
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TOTAL  RESIDUAL/I  l.OOOE-OA)  TOR  CP  IS  2.607E*10 
TOTAL  RESIOUAL/I  l.OOOE-061  FOR  HI  IS  2.63SE*10 
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TOTAL  RESIOUAL/I  l.OOOE-06)  FOR  R1  IS  3.18SE*04 
TOTAL  RESIDUAL/I  l.OOOE-OAl  FOR  Ul  IS  «.0e8E*0A 


TOTAL  RES10UAL/(  l.OOOE-OA)  FOR  VI  IS  4 

TOTAL  RESlOUAL/i  l.OOOE-04)  FOR  N1  IS  8 

TOTAL  RESIOUAL/I  1.000E>0«)  FOR  KE  IS  2 

TOTAL  RESIOUAL/(  1.000E>04I  FOR  EP  IS  2 

TOTAL  RESIOUAL/I  1.000E*04)  FOR  N1  IS  2 

NHOLE-FIELO  RESIDUALS  BEFORE  SOLUTIONS 
HHOLE-FIELO  SUH  OF  ABSIVOL.FLOH  RESIDUALS) 
HHOLE-FIELO  SUM  OF  ABSI  RESIDUALS )  OF  U1 
NHOLE-FIELO  SUM  OF  ABSI  RESIDUALS )  OF  VI 
HHOLE-FIELO  SUM  OF  ABSI  RESIDUALS)  OF  HI 
HHOLE-FIELO  SUM  OF  ABSI  RESIDUALS )  OF  KE 
HHOLE-FIELO  SUM  OF  ABSI  RESIDUALS )  OF  EP 
HHOLE-FIELO  SUM  OF  ABSI  RESIDUALS )  OF  HI 
«  SUMS  HAVE  BEEN  DIVIDED  BY  RESREFINAME) 
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AT 

PATCH 
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SOURCE 

OF 
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AT 

PATCH 

NAMED: 

SMALLS 

NET 
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AT 

PATCH 
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MULL4 
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SOURCE 

OF 
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AT 

PATCH 

NAMED: 

SHALL4 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED: 

MULLS 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAXD: 

SMALLS 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED: 

MULL4 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED: 

SMALL4 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED: 

MULL7 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED: 

SMALL7 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED: 

NNALL8 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED: 

SMALLS 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED: 

SHALL9 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED: 

HMALLXO 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED: 

UULLXn 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED: 

HMALLXX 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED: 

LMALLXX 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

fUMED: 

HMALLX2 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED; 

LMALLX2 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED: 

HMALLXS 

MET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED: 

LMALLXS 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED: 

NHALL14 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED: 

LNALLX4 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED: 

KHALLXS 

NET 

SOURCE 

OF 

UX 

AT 

PATCH 

NAMED: 

LMALLXS 

I 


I 


144E*04 

S2IE«04 

930E*08 

430E*10 

430E*I0 

S.1B5E»04 

4.088E»04 

4.144E*04 

B.S2SE04 

2.930E«^08 

I.430E+10 

2.430E*10 


*-4.494E»00 

■  1.728E-01 
•-1.27SE-02 
•-3.I72E-01 
■-2.X79E-02 
•-1.7S7E-OX 
■-2.X94E-02 
■-2.X44E-0X 
■'X.X02E-02 
•-X.940E-0X 
•-4.524E-03 
•-X.8X5E-0X 
»  2.277EtOX 
•  O.OOOEt^OO 
«  B.549E-01 

■  X.44XE*00 
•-2.399E-02 
•-9.981E-02 
S-4.374E-04 
S-X.000E400 
•-X.B9SE-02 
•-6.707E-0X 
»-2.47'c-02 
»-X.2fc4E-01 
■-X.2X0E-02 
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NET  SOURCE  OF  U1  AT  PATCH 
NET  SOURCE  OF  U1  AT  PATCH 
NET  SOURCE  OF  U1  AT  PATCH 
NET  SOURCE  OF  U1  AT  PATCH 
NET  SOURCE  OF  U1  AT  PATCH 
NET  SOURCE  OF  U1  AT  PATCH 
NET  SOURCE  OF  VI  AT  PATCH 
NET  SOURCE  OF  VI  AT  PATCH 
NET  SOURCE  OF  VI  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VI  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VI  AT  PATCH 
NET  SOURCE  OF  VI  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VI  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 
NET  SOURCE  OF  VX  AT  PATCH 


NAHEOi  HHAXLXA  ■•S.4A5E-0X 
NAHEOi  LNALLXA  ■-X.209E-02 
NAMED:  MHALLX7  B-S.424E-0X 
NAMED:  LNALLX7  ■-S.472E-03 
NAMED:  HHALLX8  b>S.S44E-0X 
NAMED:  LHALLX8  >  8.X09E-D2 
NAMED:  HHALLXO  ■  0.00OE*OO 
NAMED:  LHALLXD  ■•X.XSDE^OO 
NAMED:  MNALLXX  >  2.250E*00 
NAMED:  LHALLXX  >  9.2S9E-0S 
NAMED:  HHAUX2  «  2.302E-0X 
NAMED:  LNALLX2  ■  4.284E-04 
NAMED:  KHAt.LX3  b.9.006E-0X 
NAMED:  LNALL13  •  1.D42E-03 
NAMED:  HNALLX4  «-l.X44E»00 
NAMED:  LNALLX4  ■  1.342E'D4 
NAMED:  HNALLX5  •-2.973E*0X 
NAMED:  LHALLX5  >-4.722E-07 
NAMED:  KHALLX4  B-S.AEeE'DX 
NAMED:  U4ALLX6  •-4.430E>D7 
NAMED:  HHAULX7  •-9.907E>D1 
NAMED:  LHALLX7  »4.934E-07 
NAMED:  HHALLX8  ■-X.293E*00 
NAMED:  LNALLIS  •-X.4X5E-0X 
NAMED:  EHALLX9  •  ft.ibSE'OX 
NAMED:  HHALLX9  *  S.758E>0X 
NAMED:  ENALL20  •-e.75ZE>03 
NAMED:  NHALL20  *  3.S21E>0Z 
NAMED:  EHALLZX  *  1.254E-0X 
NAMED:  »*(ALLZX  •  1.02ZE-0X 
NAMED:  EHALLZZ  «  5.847E-01 
NAMED:  NNALLZZ  •  1.480E-0X 
NAMED:  EHALLZ3  «  1.075E-01 
NAMED:  NNALL23  «  9.808E-0Z 
NAMED:  EHALLZ4  ■  X.b94E-0X 
NAMED:  FMALLZ4  ■  9.045E-0Z 
NAMED:  EHALXZ5  «  X.Xb6E-0X 
NAMED:  fMALLZS  *  X.5X8E-DZ 
NAMED:  ENAILZ6  *-Z.047E«0X 
NAMED:  HHALLZ6  «-3.995E-0X 


NCr  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 
NET  SOURCE  OF 


NX 

AT 

PATCH 

MX 

AT 

PATCH 

NX 

AT 

PATCH 

NX 

AT 

PATCH 

MX 

AT 

PATCH 

MX 

AT 

PATCH 

MX 

AT 

PATCH 

NX 

AT 

PATCH 

NX 

AT 

PATCH 

MX 

AT 

PATCH 

MX 

AT 

PATCH 

NX 

AT 

PATCH 

NX 

AT 

PATCH 

NX 

AT 

PATCH 

MX 

AT 

PATCH 

NX 

AT 

PATCH 

NX 

AT 

PATCH 

MX 

AT 

PATCH 

NX 

AT 

PATCH 

NX 

AT 

PATCH 

NX 

AT 

PATCH 

NX 

AT 

PATCH 

NX 

AT 

PATCH 

NX 

AT 

PATCH 

NX 

AT 

PATCH 

MX 

AT 

PATCH 

NX 

AT 

PATCH 

NX 

AT 

PATCH 

NX 

AT 

PATCH 

NX 

AT 

PATCH 

RX 

AT 

PATCH 

R1 

AT 

PATCH 

RX 

AT 

PATCH 

RX 

AT 

PATCH 

RX 

AT 

PATCH 

RX 

AT 

PATCH 

RX 

AT 

PATCH 

RE 

AT 

PATCH 

RE 

AT 

PATCH 

RE 

AT 

PATCH 

NAMED:  JETOUT 
NAMED:  MULLS 
NAMED:  SNALLS 
NAMED:  MULL4 
NAMED:  SNALL4 
NAMED:  MULLS 
NAMED:  SNALLS 
NAMED:  MULL6 
NAMED:  SNALLS 
NAMED:  MULL? 
NAMED:  SHALL7 
NAMED:  MULLS 
NAMED:  SNALLS 
NAMED:  SNALL9 
NAMED:  EHALL19 
NAMED:  MULLS* 
NAMED:  EHALL20 
NAMED:  MULLED 
NAMED:  ENALL21 
NAMED:  MULL21 
NAMED:  ENALL22 
NAMED:  MULL22 
NAMED:  EHALL23 
NAMED:  MULL23 
NAMED:  ENALL24 
NAMED:  MULL24 
NAMED:  EHALL25 
NAMED:  MULL25 
NAMED:  EHALL26 
NAMED:  MULL26 
NAMED;  UETIH 
NAMED:  JETOUT 
NAMED:  OUTSIDE 
NAMED:  TOPOUT 
NAMED:  TOPP 
NAMED;  FRONT 
NAMED:  ENOENO 
NAMED:  RESOURCE 
NAMED:  JETOUT 
NAMED;  MULLS 


>  4.174E»04 
■-l.SXSE«02 
■-l.OISE^OO 

■  S.301E-03 
•-1.247E-06 

■  1.930E-03 
■-S.121E-09 
■-S.027E-03 
a-1.450E-0S 
•-S.74IE-0S 
•-X.502E-0S 
*-1.9SlE-D2 
■-1.S4DE-DS 
«-S.92SE»01 
b-2.02SE*00 
B  S.57SE-02 
B  4.467E-D2 
B-1.278E-0I 
B  2.422E>01 
B-X.i24E-0I 
B-2.73SE-01 
B-S.240E-01 
B  X.6SXE-02 
■•X.SSOE'Ol 
B  9.4X4E-02 
S-4.5S8E-02 
B  4.SS2E-04 
B  X.09SE>0X 
B-2.524E*01 
B-X.X40E-01 
b»5.849E»01 
B  3.649E*0X 
B  4.X5SE»0Z 
b-4.010E«02 
«  2.728E*0X 
B  X.859E*01 
b-4.022E«01 
B  4.X00E«05 
B  S.46SE»D2 
■-X.0e9E*10 


NET 

SOURCE 

OF 

KE 

AT 

PATCH  NAfSO:  SHALL3 

«-l.SD8E»10 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  MULL4 

■  3.850E»D9 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  SHALL4 

■-D.lllE^lO 

•NET 

SOURCE 

OF 

KE 

AT 

PATCH  NAMED:  »MALL5 

■  4.497E*09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH  NAMED:  SHALLS 

«-S.297E»lD 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  IMALLA 

>  4.514E*D9 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  SHALLS 

■•1.940E*11 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  M4ALL7 

a  l.eSSE«09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  SHALL? 

■-2.4J7E»11 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  »MALL8 

a  S.D60E«09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  SHALL8 

a-2.84SE*ll 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  SHALL9 

a-S.320E«lD 

NET 

SOURCE 

OF 

KE 

AT 

PATCH  NAMED:  HNALLID 

a  D.DODE^-DO 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  LMALLID 

a-2.D0SE«10 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  HHALLll 

a  8.344E»10 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  LNALLIl 

a-3.5S7E*08 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  HHAIL12 

a  1.310E*10 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  LHALL12 

a-1.620E*07 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  HNALL13 

a  7.155E*10 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  LNALL13 

a>8.321E«09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  HNALL14 

a  9.X99E»10 

NET 

SOURCE 

OF 

KE 

AT 

PATCH  NAMED:  LNALL14 

■-4.843Et09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH  NAMED:  HNALL15 

a  4.299E«10 

NET 

SOURCE 

OF 

KE 

AT 

PATCH  NAMED:  LNALL15 

a-4.410E«09 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  HHALLIS 

a  1.840E«^11 

NET 

SOURCE 

OF 

KE 

AT 

PATCH 

NAMED:  LNALL14 

a-2.S95E»09 

NET 

SOURCE 

OF 

KE 
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APP£3«DXX  C 
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AIRFLOW  CHARACTEHIZATIOM  SHEFT 


meeilAllt  T««CiP8  Ot  tUyy  Aircraft  In  Archwd  Structure  Hu«h  House 
LOCATlONt  lUXly  A.F.B 
OATSl  aa 

AIKCNArTi  F»4H 

AnaiME.  J-7<>aCE-a  **»•»•  »tF«i»*lOI0  (IB.  ■B.i; 

EVENTi  5 

CONDITIONSt  Power.  Starboard  Engine  -  ldle_P.ewfer - 


AUOMENTEB  IMkET 


AIRFLOW  CHARACTERIZATION  SHEET 


Test  PROCRAM:  Thesis 

kocATioNi  Naval  Air  Station.  JacksonviTIp 

OATCt 

AIRCRAPTi , 


Flnrit^a 


'n/a' 


gwoiMEi  J-79-GE-8 _  aiPRKtsiOM  (la.  ara.i;  M/a 

CVCNTi  _ 


CONDITIONS:  Single  Engine  with  A/B  g  centerlinp 


APP£34DXX  J> 


EXHAUST  FLOW  DATA  SHEET 


vaav  BBAABAu.  TMCing  of  lUvy  Alrcr«fc  In  Archnd  Scruccur*  Hush  House 
irCATlOH?  luily  A.F.B.  f  ilkn  Mcnnlo.  tnxns  ~ 

DATE:  9  Moveaber  1982 

AIRCRAFT:  _ 

RWCIWRI  J-79-CE-8 

KVENT:  S _ 

CONBITIONS:  Port  Engine  -  A/B.  SesrboArd  Engine  -  Idle  Power 


AUCMCNTCR  TURE 


Tavg=  454 

augmentation  AATIOS  S.3 


L28 


EXHAUST  FLOW  DATA  SHEET 


re»r  prooiiam:  Thesis _ ■ 

LOCATION:  Naval  Air  Station,  Jacksonville.  Florida 

OATE: 

AINCNAFTi 

ENGINE:  J-79-GE-8 
EVENT:  _ 

coNoiTiONS:  Single  Engine  with  A/B  g  centerline 


l’ 


I 


I 


AUOMENTEN  TUNE 


FLOW  IINak/***.)*  2373 
AUOMCNTATION  NATiO>  iQ  .  6 


Vavg=  235.5 
Tavg=  401.6 


I 
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APPE34DXX  E 

l>A.*rA  TRANSFER  BV  TAPE 


The  IBM  3033  main-frame  presently  does  not  have  "PHOTON" 
graphics  capabilities.  In  order  to  produce  various  visual 
plots  that  enhance  the  interpretation  of  tabulated  data,  the 
"DF09"  and  "DF12"  files  must  be  properly  formatted  and  then 
transferred  to  tape  (for  portability)  allowing  use  of  "PHOTON" 
on  another  computer  which  is  "PHOTON"  capable. 

The  computer  center  will  provide  a  non-labelled  tape  for 
data  transfer.  The  operator,  when  queried  by  the  programs 
below,  will  mount  the  tape  allowing  the  "DF09"  and  "DF12" 
files  to  be  transferred  to  the  tape. 

To  prepare  the  files  for  transfer  to  tape  a  copy  is  made 
to  ensure  that  the  original  data  file  is  safe.  Program  "FORM09 
JCL"  properly  formats  and  then  writes  the  "DF09"  file  to  a  new 
file  called  "DF09C".  Program  "FORM12  JCL”  properly  formats  and 
then  writes  the  "DF12"  file  to  a  new  file  called  "DF12C".  The 
data  sets  "DF09C"  and  "DF12C"  which  are  produced  using  MVSHELP 
must  be  at  least  as  large  as  the  original  "DF09"  and  "DF12" 
data  sets  and  created  before  running  "FORM12  JCL"  or  "FORM09 
JCL"  . 

With  the  two  files  now  properly  formatted  the  file  "DSTP 
JCL"  is  used  to  transfer  the  designated  file  to  tape.  In  the 
following  example  the  "DF09C"  file  is  transferred  as  the 


130 


second  file  to  the  tape  called  "NICOLA".  File  "TSCAN  JCL"  then 
scans  the  tape  to  ensure  the  program  was  actually  transferred 
to  the  tape. 

Two  other  files,  "TAP  JCL"  and  "DUMP  JCL",  can  dump  a 
file  from  tape  to  the  virtual  machine  and  return  the  first 
block  of  the  tape  file  respectively.  With  the  files  on  tape 
and  checked  the  tape  can  be  dismounted  by  the  operator  at  the 
computer  center  and  transported  to  a  "PHOTON"  capable  system 
in  the  "formatted"  form. 

The  various  files  referenced  are  now  listed. 


FILE:  OSTP  JCL  A1 


//OSTP  JOB  <1541, 9999),' ERIC SCLASSsE 
//  EXEC  PGn=IEBGENER 
//SYSPRIMT  00  SYSOUT=A 
//SYSIN  00  OUMMY 

//SYSUTl  00  0ISP=SHR,0SN:nsS.S1541.PHOENICS.0FO9C 
//SYSUTZ  00  UNIT=3<tOO-6,VOLsSERsNICOLA,DISP=INEH,PASS), 

//  LABEL=(2,NL,,0UT), 

//  DCB:(RECFH=FB,LRECL=B0,BLKSIZEsB00,0EN=S,OPTC0=Q) 


I 


I 
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//TAP  JOB  1 1541 >9999 )> ‘TAP  JCL‘>CLASS>E 
//•HAIN  CAROS«(SO) 

//  EXEC  PGM^ZEBGENER 
//SYSPRINT  DO  SYSOUT«A 
//SYStN  DO  DUMMY 

//SYSUTl  OD  UNITsS400-5,VOL«SER>NICOLA>OISP«l OLD, PASS )> 
//  LABEL«(2>NL>>IN)> 

//  OCB>( RECFMxFB>LRECL>80>BLKSIZE«800,OPTCDbB  I 

//SYSt/TZ  DO  SYS0Ur<B>0C8«BUCSZZE*80 
// 


FILE:  DUMP  JCL  A1 


//DUMP  JOB  <1541, 9999)> ‘DUMP*, CLASS=E 
//  EXEC  TAPE.VOUNsNICOLA 
//SYSIN  00  « 

OMPREC< 10,272,1)  ,SKPFIL« 10,1)  ,XECUTE ( 999 ) 
/• 

// 


FILE:  TSCAN  JCL  AI 


//6SCAN  JOB  1 1541, 9999), ‘NICOLA’, CLASSsE 
//  EXEC  TSCAN,VOLIN=NICOLA 
// 


FILE;  F0ltM09  JCL  A1 


//F0«tM09  JOB  ( 1541,9999  )> ‘CHAM*  >CLASB>B 
//«MAIN  LINES«(20) 

//  EXEC  FORTVCLG 
//FORT.SYSIN  OD  • 

C  PROGRAM  F0RM09 
C 

C  Incroasa  PARAMETER  ••L'tinBs  if  aor*  variablM  or  •  largar  grid 
C  is  rsquircd. 

C 

PARAMETER( NPHIP*eO ) 

C  MB  HKYMAX  susi  b*  >*  NZ  and  l«<aNy 
PARAMETER! NXYMAXcZOOOO ) 

CHARACTER«4  MAMEINPHIP),MESS(10) 

CHARACTERaiS  FILIN, FILOUT 
CHARACTER*!  ANS 
CHARACTERalO  ACCESS 
LOGICAL  LBUF(10),ST0RE(NPHIP),DEBUG 
DIMENSION  RBUF(NXYMAX),1BUF(10) 

C 

100  FORMAT! IX, 71 10) 

101  FORMAT! 61 1PE13.6)) 

102  FORMAT! IX, 79L1) 

103  FORMAT!  IX,  19A4) 

C 

C  Rasei  DEBUG*. TRUE,  io  gsi  arror  aassagas.. 

DEBUG*. TRUE. 

C 

C  Sat  Logical  Units... 

LUINa9 

LUOUT=10 

LUBUG:6 

C 

C  Set  File  Names... 

C  Name  of  unformatted  restart  file 
C  HRITE!6,*)  •  RESTART  FILENAME?’ 

C  READ!5,'!A)')  FILIN 

C  Name  of  formatted  restart  file 
C  MRITE!6,*)  •  FORMATTED  FILENAME?’ 

C  READI5,'!A)' )  FILOUT 

C 

C  Direct  access  or  scc^Jential? 

C  IHRITE!6,«)  *  DIRECT  ACCESS  IPHIDA)  FILE?  Y/N' 

C  READ!5f ‘!A1)' )  ANS 

C  IF! ANS. EG. 'Y'. OR. ANS. EG. ’Y’)  THEN 

C  Version  1.4  default 
C  ACCESS* 'DIRECT' 

C  LENREC*256 

C  ELSE 

C  Version  1.3 

ACCESS* 'SEGUENTIAL' 


I 


I 


! 
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oo  no  n  o  o  o  o  o  n  f>  o 


FILE:  FORM09  JCL  A1 


LENREC^lOO 

ENOIF 

HRITEI6i«)  ’  DEFAULT  LENREC?  flOO  VI. 3,  ZS6  VI. A)  -  Y/N* 
READIS,MA)' )  ANS 

IF(ANS.NE.'Y'.ANO.ANS.NE.*Y')  THEN 
MRITEI6,*)  ‘  LENREC?  CHORDS  NOT  BYTESI* 

REA0IS,‘(1)’ )  LENREC 
ENDIF 

NBYTES  i«  ‘th*  r«Mb«r  of  byfM  por  singlo  procisien  iwrd. 

NBYTESS4 

Rocord  length  in  bytes. 

LREC sNBYTES«LENREC 

OPEN  input  file: 

OPEN! UN1T-9.F0RH« ‘UNFORMATTED* .STATUS- ‘OLD ' > 

1  ACCESSsACCESS.I0STAT-IERR,ERR*90) 

IFIIERR.NE.O)  GO  TO  90 
IFCACCESS.EQ. ‘DIRECT* )  THEN 
READCLUIN.RECsl)  IMESSC I  ),I-I>10) 

ELSE 

REAOCLUINI  (KESS(II>I-1,10) 

ENDIF 

IFC DEBUG)  MRITE(LUBUG.«)  HESS 
IFCACCESS.EQ. ‘DIRECT*  )  THEN 
READtLUlN.REC«2)  I LBUFC I ).I«I>4) 

ELSE 

READ(LUIN)  (LBUF(I)>I-1,4) 

ENDIF 

IFCOEBUG)  NRITE(LUBUG>«)  (LBUFCI  )>I=I,4) 

IFCACCESS.EQ. ‘DIRECT*  )  THEN 

READCLUIN,REC=S)  NX.NY, NZ.NPHI. I IBUFC I ), 1-1.6 ) .LENREC 
ELSE 

READCLUIN)  NX.NY.NZ.NPHI.CIBUFCI).I*1.6). LENREC 
ENOIF 

NXNY=NX«NY 

IFCNXNY.GT.NXYMAX)  GOTO  92 
IFCNPHI.GT.NPHIP)  GO  TO  93 
IFC  DEBUG)  NRITECLUBUG.*)  tCC.NY.NZ 

OPEN  output  file:  FORMATTED.  SEQUENTIAL  end  of  FIXED  RECORD  LENGTH. 
OPENC UNIT-10, FORM- ‘FORMATTED* .STATUS- ‘ OLD  * . 

1  ACCESS- ‘SEQUENTIAL ‘ ,I0STAT-IERR,ERR-9l ) 

IFCIERR.NE.O)  GO  TO  91 
HRITECLUOUT,103)  (MESSC I ).I-1.10 ) 

HRITE(LU0CJT,102)  C  LBUF C I ) .1-1 .4  ) 

NRITE C  LUOUT , 100 )  NX ,NY ,NZ ,NPHI . C IBUFI I )  ,I -1 ,6  ) , LENREC 
IFCACCESS.EQ. ‘DIRECT* )  THEN 
READCLUIN,RECs4)  RlhtCER 
ELSE 

READCLUIN)  RINNER 


FILE :  FORMO«*  JCL 


A1 


ENDIF 

NRITE(LU0UT,101)  RINNER 
IFI ACCESS. EQ. 'DIRECT' )  THEN 
READILU1N>REC«5)  (NAHE(X),I«1>NPH1) 
ELSE 

REAO(LUIN)  «NAHEIX),1.1,NPHI) 

ENDIF 


C  Sci  C36...  iwMS  for  variablos  SI  fo  NPHI  for  vorsion  1.4 
C  control  charactor*  otel  (Unloss  olroody  stored) 

1 F  ( ACCESS .  E(J . '  DIRECT ' .  AND .  NPNI  .6T .  50 )  THEN 
REA0(LUIN.REC>10I  (ST0RE(I|,Is1,NPHII 
00  6  I>51>NPHI 


» 


J=I-1S 


1F( .NOT.STORE(I))  THEN 

if(i.lt.ioo)  then 

HRITE(NAriE(I),'(  "C"»I2,"  ")')J 
ELSE 

NRITE(NAME(I),'(  "C"»I5)' )  J 
ENDIF 
ENDIF 

A  CONTINUE 
ENDIF 


WRITE) LUOUT, 105)  (NAMECI ),1«1>NPHI ) 

IF(OEBUG)  WRITE! Ll»UG,o)  (NAME) X ),I*1,NPHI ) 
IF(ACCESS.Eq. 'DIRECT' )  THEN 
READ(LUIN,REC«6)  (RBUF(I),I«1,N>() 

ELSE 

READ(LUIN)  )RBUF(1),I.1,NX) 

ENDIF 


WRITE(LU0UT,101)  (RBUF(I  ),Isl>»«<) 
IF ( ACCESS. EQ. 'DIRECT')  THEN 
READ(LUIN,REC=7)  (RBUF(I),m>NY( 
ELSE 

REAO(LUIN)  (RBUF(I),Isl,NY) 

ENDIF 

WRITE)  LUOUT, 101)  ( RBUF)  I  ),in»NY ) 
IF) ACCESS. EQ. 'DIRECT' )  THEN 
READ)LUIN,RECs8)  (RBUF)I),Is1,NZ) 


REAO(LUIN)  (RBUF)I),Is1,nz) 

ENDIF 

WRITE) LUOUT, 101)  (RBUF)I ),Icl,NZ) 
IF)ACCESS.Eq. 'DIRECT' )  THEN 
REA0)LUIN,RECa9)  ( RBUF(  I ),I«1,NZ) 


READ)LUIN)  (RBUF)I),Is1,NZ) 

ENDIF 

WRITE)  LUOUT,  101)  (RBUFd  ),m,NZ) 

IF) ACCESS. EQ. 'DIRECT' )  THEN 

READ) LUIN,REC*10 )  ( STORE) I ),I:1,NPHI ) 
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FILE:  F0RH09  JCL 


A1 


ELSE 

REAOILUZN)  (STORE(I)>I>l>NPHII 
ENDIF 

MRITEILUOUT.IOZ)  (STORE! Z )»I>1>NPHZ ) 

ZF(OEBUS)  HRITE(LUBUS,102)  (STORE! Z ),Z>1,NPHI ) 

NVARsO 

DO  5  I>Z,NI>HI 
IFISTOREd))  NVAR«NVAR«1 
S  CONTINUE 
ZREC>10 
DO  10  K*1>NZ 
DO  20  I<1>NVAR 
I FI  ACCESS. EQ. ‘DIRECT*)  THEN 
ZREC>ZREC«1 
NVREC=N»(Y/LENREC 

IF(l100(N>a(V,LENREC).NE.0)  NVREC<NVREC*1 
DO  30  IVREC>1,NVREC 
JF3(IVREC-1)«LENREC  *  1 
JL«IVREC»LENREC 
JL>HINOIJL,NXNY) 

READ! LUIN>REC«IREC )  ( RBUFI J ), J«JF , JL ) 

SO  CONTINUE 
ELSE 

READ(LUZN)  IRBUF(U)>J>1>NXNY) 

ENDIF 

HRITEILUOUT.lOl)  (RBUFI J),J«1>NXNY) 

20  CONTINUE 
10  CONTINUE 

IF( DEBUG)  MRITE(LUBU6>1900) 

1900  FORMAT I IX >■  NORMAL  STOP  IN  PROGRAM’) 

GO  TO  999 

C  Error  trapping. , . 

90  IFI DEBUG)  HRITEI LUBUG>1901 ) 

1901  FORMAT!  IX, '  ERROR  IN  OPENING  INPin*  FILE '  ) 

GO  TO  999 

91  IFI DEBUG)  HRITEI LUBUG, 1911) 

1911  FORMAT! IX, 'ERROR  IN  OPENING  OUTPUT  FILE') 

GO  TO  999 

92  IFINXNY.6T.)«<YNAX)  HRITEI  LUBUG, 1921 ) 

IF(NZ.GT.NXYHAX)  HRITEI LUBUG, 1925 ) 

1921  FORMAT! IX, 'SETTING  OF  NX»Ny  TOO  LARGE)  RESET  IN  FORM09  AND 
1,/, IX, 'RE-CREATE  THE  LOAD  MODULE') 

1923  FORMAT! IX, 'SETTING  OF  NZ  TOO  LARGE)  RESET  IN  FORM09  AND  ' 
1,/, IX, 'RE-CREATE  THE  LOAD  MODULE') 

60  TO  999 

93  HRITEI LUBUG, 193) 

193  FORMAT! IX, 'SETTING  OF  NPHI  TOO  LARGE)  RESET  IN  FORM09  AND  ' 
1,/, IX, 'RE-CREATE  THE  LOAD  MODULE') 

999  CONTINUE 
C  CLOSE  filas  ... 


FILE:  F0KM12  JCL 


A1 


//F0RM12  JOB  ( 1541 > 9999 ),*CHAH‘>CLASS>B 
//•MAIN  LINES* (20) 

//  EXEC  FORTVCLG 
//FORT  SYSIN  OD  • 

C  Fil«  n«iie:  F0RH12  lact  — ndad  15  October  1906 
C  PROGRAM  F0RM12 

C - Program  to  roHrita  a  0F12  file  into  a  formatted  READCO  file. 

C  (JPE  t  JFO  ) 

C 

C  Reset  NXYIP  if  larger  grid  required. . 

PARAMETER( NXY1P*10000 ) 

DIMENSION  AREC( NXYIP) 

C  Reset  fixed  record  length  of  0F12  file  if  required. . . . 

CHARACTER«30  FILIN, FILOUT 
LOGICAL  DEBUG 
INTEGER  SECTOR 
DATA  LENREC/100/ 

C 

100  FORMAT! 315) 

101  F0RMAT(5(1PE13.6)) 

C 

C  SECTOR  represents  the  largest  record  length  that  can  be  OPENed. 
SECTOR=6SS36 

C  KBYTES  is  the  number  of  bytes  per  single  precision  word. 

)BYTES*4 

LRECsMINOC  NBYTES«LENREC .SECTOR-MBYTES ) 

IBLK*LREC«NBYTES 

C 

C  Set  Logical  Uhits... 

LUIN*12 

LUOUT=13 

LUBUG*6 

C  Reset  DEBUG  if  error  messages  required. . 

DEBUG*. TRUE. 

C  Set  file  rmimes 

C  FILIN='0F12  • 

C  FILOUT='DF12FORM  ’ 

C 

C  OPEN  ir^t  file:  DIRECT  ACCESS,  UNFORMATTED  and  of  FIXED 
C  RECORD  LENGTH. 

C 

OPEN! UNIT*12,FORN* 'UNFORMATTED ' .STATUS* 'OLD ' , 

1  ACCESS* 'DIRECT ' ,RECL*LREC , 

1  I0STAT*IERR,ERR*90) 

IF(IERR.NE.O)  GO  TO  90 
C 

READILUIN,REC*1)  KDCPl.NYPl.NZPl 
NXY1*NXP1*NYP1 

IF!NXY1.GT. NXYIP)  GO  TO  92 
NREC*NXY1/LENREC 


non 


FILE:  F0mi2  JCL 


i 


A1 


ZSHORT-W(Vl-LENREC«NREC 

IFIISHORT.CT.O)  NREC«NRECn 
IFIISHORT.EQ.OI  ZSHOIiT«LENREC 

ZF(DEBUG)  MRITE(LUBU6»«)  NXPl>NYPl»NZPl>»e(Yl,NREC,ISHORT 

OPEN  output  filo:  FORMATTED.  SEQUENTIAL  and  of  FIXED  RECORD  LENGTH. 

OPENI  UNIT*13  .FORM=  ‘  FORMATTED  *  .STATUS*  'OLD ' . 

I  ACCESS* 'SEQUENTIAL* .IOSTAT*IERR.ERR«91 ) 

IF(IERR.NE.O)  GO  TO  91 
HRITE( LUOUT.lOO )NXP1.NTP1.NZP1 
NR*1 

DO  5  K*1,NZP1 
DO  6  1C*1.S 
IHAX*LENREC 
DO  7  JR*1,MIEC 
ZFIJR.EQ.NREC)  ZMAX*1SH0RT 
NRsMttI 

»ffiA*IJR-l )»LENREC 

READ(LUIN>REC*NR)  (AREC(NRA*1 1.1*1. ZMAX) 

7  CONTINUE 
C  Hrito  data  to  output  filo 

MRITE(  LUOUr.lOl )( ARECI II  ).II*1>NXY1 ) 

6  COriTINUE 
5  CONTINUE 

IF(OEBUG)  MRITE(LUBU6>1900) 

1900  FORNAT(  IX, 'NORMAL  STOP  IN  PROGRAM') 

GO  TO  999 

C  Error  trapping... 

90  IF( DEBUG)  KRITEf LUBUG.1901 ) 

1901  FORMAT)  IX, 'ERROR  IN  OPENING  INPUT  FILE') 

GO  TO  999 

91  IF( DEBUG)  WRITE) LUBUG, 1911) 

1911  FORMAT) IX, 'ERROR  IN  OPENING  OUTPUT  FILE') 

GO  TO  999 

92  IF) DEBUG)  WRITE) LUBUG, 1921 ) 

1921  FORMAT) IX, 'GRID  TOO  LARGE j  RESET  NXYIP  AND  RE-CREATE', 

1  /, IX, 'LOAD  MODULE') 

999  CONTINUE 
C  CLOSE  filas  ... 

CLOSE ( LUIN, STATUS* 'KEEP '  ) 

CLOSE ) LUOUT , STATUS* 'KEEP ' ) 

STOP 

END 

//GO.FT12F001  DO  DSN*HSS.S1541.PH0ENICS.DF12.DISP*)0LD ) 

//GO.FT13F001  DO  DSN=HSS.S1S41.PH0ENICS.DF12C,DISP*)0LD) 

/* 

// 
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